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ABSTRACT 


The  construction  of  an  apparatus  for  the  condensation 
performance  testing  of  a  horizontal  bundle  of  four  tubes  with 
various  refrigerants  was  completed.  The  apparatus  was 
instrumented,  and  data  reduction  software  was  developed  to 
provide  bundle  and  single  tube  condensation  data. 

Two  tube  bundles  were  tested,  smooth  copper  tubes  and  low 
integral-fin  copper-nickel  tubes,  with  two  refrigerants,  R-114 
and  R-113.  An  enhancement  ratio  of  about  2.0  for  the  overall 
heat  transfer  coefficient  was  demonstrated  for  the  finned 
tubes  over  the  smooth  tubes. 

Internal  contamination,  possibly  due  to  a  breakdown  of  the 
refrigerant  molecules  when  subjected  to  high  temperatures  in 
the  boiling  chamber,  inhibited  further  meaningful  data 
collection.  Recommendations  for  improvement  of  the  test 
apparatus  are  made. 
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NOMENCLATURE 


Aef  Effective  outside  area  of  the  tube  (m2) 

Af  Actual  area  of  a  finned  tube  (m2) 

Ai  Tube  inside  area  (m2) 

A0  Tube  outside  area  (m2) 

Ar  Surface  area  of  tube  at  the  base  of  the  fin  (m2) 

Ci  Inside  correlation  coefficient 

0o  Outside  correlation  coefficient 

Cp  Specific  heat  of  coolant  (J/kg  K) 

Deq  Equivalent  diameter  of  finned  tube  (m) 

D*  Inside  diameter  of  test  tube  (m) 

D0  Outside  diameter  of  test  tube  (m) 

Dr  Root  diameter  of  finned  tube  (m) 

g  Acceleration  due  to  gravity  (m/s2) 

L  Length  of  the  condenser  tube  (m) 

LMTD  Log  mean  temperature  difference  (K) 

N  Number  of  the  tube  from  top  of  bundle 

Ph  Phase  Change  Number 

Pr  Prandtl  Number 

Q  Heat-transfer  rate  (W) 

Q"  Heat-flux  (W/m2) 

Re  Reynolds  Number 

Rm  Tube  wall  thermal  resistance  (m2  K/W) 

Res  Swirl  Reynolds  Number 
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T^  Coolant  inlet  temperature  (K) 

Tc0  Coolant  exit  temperature  (K) 

Tsat  Vapor  saturation  temperature  (K) 

U0  Overall  heat-transfer  coefficient  (W/m2'K) 

Greek 

Inside  heat-transfer  coefficient  (W/m2  K) 

at  Average  inside  heat-transfer  coefficient  (W/m2  K) 

a0  Outside  heat-transfer  coefficient  (W/m2  K) 

aQ  Average  outside  heat-transfer  coefficient  (W/m2  K) 

Ahv  Specific  enthalpy  of  vaporization  (J/kg) 

AT  Temperature  rise  of  coolant  across  condensing 

length  (K) 

A  Thermal  conductivity  (W/m  K) 

r  Mass  flow  rate  of  coolant  (kg/s) 

p  Density  (kg/m3) 

ef f ~f  Fin  efficiency 

r)  Dynamic  viscosity  (N  s/m) 


I. 


INTRODUCTION 


A.  BACKGROUND 

With  the  advent  of  more  complex  shipboard  weapons/combat 
systems  and  increasing  fuel  costs,  the  United  States  Navy  has 
increasingly  recognized  a  need  for  an  energy-efficient,  light¬ 
weight,  and  high-capacity  air-conditioning  system.  Such  an 
advanced  system  has  been  proposed  and  is  currently  in  its  test 
and  evaluation  phase  for  the  DDG-51  program  [Ref.  1] . 

Among  the  substantial  differences  in  the  advanced  air- 
conditioning  plant  over  extant  systems  is  the  use  of  titanium 
finned  tubes  in  the  condenser,  in  place  of  the  smooth  copper- 
nickel  alloyed  tubes  previously  used.  The  use  of  titanium  in 
the  air-conditioner's  condenser  presents  two  significant 
advantages.  First,  since  the  United  States  Navy  uses  ambient 
seawater  as  a  heat  sink,  the  corrosion  resistance  of  titanium 
should  measurably  improve  system  integrity  and  reliability. 
Secondly,  the  advanced  air-conditioning  system  expects  to 
realize  a  weight  savings  per  unit  of  2286  kg  (5040  lbs)  and 
the  major  portion  of  this  weight  savings  will  come  from  the 
use  of  the  lighter-weight  titanium  condenser  tubes  [Ref.  1]. 

Since  the  early  1970 's  the  United  States  Navy  has  used  R- 
114  as  its  primary  refrigerant.  The  advantages  of  R-114  over 
the  more-widely-used  commercial  refrigerants  in  a  naval 
application  derive  from  its  low  toxicity,  temperature 
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stability,  stability  when  in  contact  with  moisture,  and  its 
applicability  to  lower-pressure  systems  [Ref.  1].  R-114.1  will 
be  the  refrigerant  utilized  in  the  advanced  air-conditioning 
plant. 

As  no  substantial  data  base  on  the  heat-transfer  capacity 
exists  for  a  system  utilizing  titanium  condenser  tube  bundles 
configured  as  proposed  in  the  advanced  air  conditioning  plant 
and  consequently,  no  way  exists  to  measure  this  design's 
performance  against  competing  designs,  the  need  arose  to 
develop  a  test  apparatus  to  provide  this  data  base  in  order  to 
accurately  predict  future  system  performance.  The  design  and 
subsequent  construction  of  this  test  apparatus  at  the  Naval 
Postgraduate  School  began  in  1987  and  is  documented  in  the 
Theses  of  LCDR  David  S.  Zebrowski  [Ref.  2]  and  of  LCDR  Thomas 
J.  Murphy  [Ref.  3],  The  test  apparatus  was  also  intended  to 
serve  as  a  test  platform  for  advanced  boiling  surfaces 
proposed  for  use  in  the  advanced  air-conditioning  plant.  The 
incorporation  of  the  advanced  boiling  surface  test  platform 
and  a  condenser  tube  test  platform  into  a  single  test 
apparatus  allows  the  widest  latitude  in  examining  the  effects 
of  various  heat  inputs,  coolant  flow  rates,  and  various  levels 
of  refrigerant  contaminants  on  individual  components  and 
overall  system  performance.  Zebrowski  and  Murphy  constructed 


JNew  refrigerants  may  have  to  be  used  in  the  future  due  to 
the  "ozone  problem." 
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the  basic  apparatus,  and  Murphy  operated  the  system  for 
preliminary  evaporator  measurements. 


B.  OBJECTIVES 

The  major  objectives  of  this  thesis  were: 

1.  Refine  the  design  and  complete  the  fabrication  of  the 
apparatus  for  the  testing  of  condenser  tube  bundles 
utilizing  R-114  as  the  working  fluid. 

2.  Develop  and  instrument  a  data-acquisition  system  with 
the  associated  software  to  provide  a  data  base  on  the 
heat-transfer  capacity  of  the  condenser  section  of  the 
test  apparatus. 

3.  Validate  the  test  apparatus  by  comparing  condenser 
performance  against  extant  data  bases  derived  from 
conventional  condenser  tube  bundles. 
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II.  LITERATURE  SURVEY 

A.  GENERAL  OBSERVATIONS 

Condensation  is  the  phase  transformation  process  in  which 
a  vapor  is  transformed  to  liquid  by  removal  of  latent  heat. 
The  promotion  of  condensation  in  heat  exchangers  is  used 
extensively  in  applications  for  propulsion  engineering  and 
air-conditioning/refrigeration  cycles.  Due  to  its  importance 
in  the  aforementioned  fields,  considerable  research  has  been 
directed  at  the  factors  influencing  the  process  and  ways  in 
which  the  process  can  be  enhanced.  Such  factors  as  various 
modes  of  condensation,  surface  orientation  to  vapor  flow,  the 
shear  forces  exerted  on  condensate  film,  various  external  and 
internal  surface  enhancements,  turbulent  effects  due  to 
cascading  condensate  flow  from  another  surface,  and  the 
effects  of  fluid  properties  on  the  process  are  detailed 
extensively  in  various  reviews.  This  thesis  deals  exclusively 
with  film  condensation  of  refrigerants  on  small  horizontal 
tube  bundles  in  a  quiescent  vapor  where  tubes  with  external 
fins  are  compared  to  smooth  tubes. 

B.  SINGLE  TUBE  INVESTIGATIONS 

1 .  Smooth  Horizontal  Tube  Studies 

The  first  comprehensive  condensation  model  was 
developed  by  Nusselt  in  1916  [Ref.  4]  based  on  the  assumption 
that  a  quiescent  vapor  at  saturation  temperature  coming  into 
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contact  with  a  wall  surface  below  saturation  temperature  would 
condense  and  form  a  continuous  film  of  condensate  growing  in 
thickness  as  the  film  flowed  off  the  surface  under  the 
influence  of  gravity.  The  condensate  film  at  the  vapor-liquid 
interface  would  be  at  the  vapor  saturation  temperature  with  a 
temperature  gradient  in  the  film  down  to  the  wall  surface 
temperature.  No  radiation  or  convection  would  take  place  at 
the  liquid-vapor  interface  as  both  liquid  and  vapor  are  at  the 
same  temperature,  however  the  amount  of  vapor  condensing 
corresponds  to  the  quantity  of  heat  flowing  through  the  film 
by  thermal  conduction.  In  the  case  of  a  horizontal  tube,  he 
assumed  that  condensate  flows  in  a  laminar  manner  around  the 
sides  of  the  tube  and  off  the  tube  in  a  continuous  sheet.  The 
heat  transfer  coefficient  would  be  maximum  at  the  top  center 
of  the  tube  decreasing  gradually  around  the  surface  of  the 
tube,  as  the  inherent  thermal  resistance  in  the  film  grows 
with  the  thickness  of  the  film,  and  eventually  goes  to  zero  at 
tube  bottom.  Nusselt's  model  is  limited  by  disregarding 
surface  tension  forces  that  tend  to  hold  up  condensate  at  the 
tube  bottom  until  overcome  by  gravity,  resulting  in  the 
production  of  condensate  droplets  rather  than  a  continuous 
sheet.  Nusselt's  expression  for  the  average  heat  transfer 
coefficient  for  a  single  tube  subjected  to  a  constant  heat 
flux  is  given  by: 

aD  =  .655[A^p^g-Ahv/nf*Do*Q"]1/3 
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Although  constrained  by  the  aforementioned  limitations, 
Nusselt's  model  remains  the  conservative  bench-mark  against 
which  all  other  models  are  compared. 

2 .  Exterior  Surface  Enhanced  Tubes 

Research  into  ways  to  improve  condensation  performance 
in  condensers  is  motivated  by  the  idea  that  any  force  that 
acts  to  thin  the  condensate  film  promotes  an  enhancement  of 
heat  transfer  by  minimizing  the  resistance  to  heat  flow. 
General  enhancement  techniques  include  the  fabrication  of 
surfaces  that  promote  dropwise  condensation,  wrapping  exterior 
surfaces  with  wire,  installation  of  porous  drainage  strips, 
and  the  fabrication  of  finned  tubes  with  fins  of  various 
geometries  and  various  spacings  along  the  tube.  This  thesis 
deals  with  a  low  integral-fin  tube  and  its  heat  transfer 
enhancement  over  a  smooth  tube. 

In  considering  a  low  integral-fin  tube  during 
condensation,  there  exists  two  distinguishable  regions  on  the 
circumference  of  the  tube;  a  flooded  region  and  an  un-flooded 
region.  The  flooded  portion  of  the  tube  defines  the 
condensate  retention  angle  of  the  tube  with  respect  to  the 
tube  circumference.  The  smaller  the  condensate  retention 
angle  of  the  tube,  the  larger  the  heat  transfer  capacity  of 
the  tube.  The  fin  can  be  divided  into  three  regions;  the  fin 
t  os,  the  sides  of  the  fins,  and  the  fin  root  area.  The 
majority  of  heat  transfer  takes  place  at  the  fin  tips. 
Surface  tension  forces  pull  the  condensate  from  the  tips  down 
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the  fin  sides  into  the  flooded  root  area,  where  gravity  drains 
the  condensate  to  the  bottom  of  Lhe  tube.  The  amount  of 
condensate  retained  along  the  circumference  of  the  tube  is 
dependent  upon  the  ratio  of  surface  tension  forces  to  gravity 
forces. 

Beatty  and  Katz  [Ref.  5],  in  1948,  developed  a 
comprehensive  model  for  the  prediction  of  the  heat  transfer 
coefficient  of  finned  tubes  based  on  their  experimental 
results  for  various  test  fluids  (including  R-22)  and  finned 
tubes  of  various  metallic  compositions,  and  various  fin 
geometries.  Their  model  assumes  gravity-dominated  flow  and 
neglects  surface  tension  effects  completely.  Subsequent 
experimental  results  from  other  sources  indicate  that  the 
Beatty  and  Katz  model  overpredicts  the  heat  transfer 
coefficient  as  surface  tension  increases  or  as  fin  density 
increases.  Nevertheless,  the  Beatty  and  Katz  model  conformed 
to  their  reported  overall  heat  transfer  enhancement  of  up  to 
2.3  for  R-22  on  finned  tubes  compared  to  smooth  tube  data.  The 
Beatty  and  Katz  model  is  given  by: 


o  o  1/4  1  -i  1/4 

ao  =  .689.[>|-0f-g-/Vn£-ATvf]  -to—l 

eq 

where. 


[J_]1/4  = 

lD  j 
eq 


ef  d; 


1/4  +1-30’  [neff-f]  [A 


Af  1 

1  biTT1 

ef  L7 


(2.2) 


(2.3) 
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and, 


L  =  tt.  (D^  _  )  /  4  •  Dq 

A  three-region  theoretical  model,  based  on  their  experimental 
results  with  three  rectangular  finned  tubes,  was  presented  by 
Karkhu  and  Borovkov  [Ref.  6],  in  1971,  including  surface 
tension  forces.  Measured  heat-transfer  coefficients 
demonstrated  a  50  to  100  percent  increase  in  vapor-side  heat 
transfer  coefficients  for  steam  and  R-113  condensing  on  finned 
tubes  compared  to  smooth  tubes.  Unfortunately,  they  did  not 
report  enhancements  separately  for  the  two  fluids.  In  1980, 
results  reported  by  Carnavos  [Ref.  7]  condensing  refrigerants 
on  various  finned  tubes  demonstrated  an  enhancement  of  up  to 
400  percent  in  the  heat  transfer  coefficient  over  smooth  tube 
results.  Work  done  by  Sauer  and  Williams  [Ref.  8],  in  1982, 
on  the  condensing  performance  of  finned  tubes  with  oil- 
contaminated  R-113  demonstrated  a  serious  degradation  of 
performance  when  the  surface  tension  to  density  ratio  was 
large.  The  conclusion  formed  was  that  the  higher-surface- 
tension  oil  remained  in  the  fin  gaps  rendering  the  finned 
surface  ineffective.  Results  reported  by  Honda  et  al.  [Ref. 
9],  in  1983,  condensing  R-113  on  finned  tubes  of  various 
geometries  showed  improvements  of  900  percent  for  the  vapor- 
side  heat  transfer  coefficients  at  a  constant  vapor-to-tube- 
wall  temperature  difference.  The  results  reported  by  Kabov 
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[Ref.  10]  in  1984  with  refrigerants  R-12  and  R-21,  indicated 
that  the  bulk  of  the  latent  heat  was  in  fact  removed  on  the 
lateral  fin  surfaces  and  that  an  optimum  fin  height  and 
spacing  was  dependent  upon  the  ratio  of  surface  tension  to 
gravity  forces.  The  results  published  by  Masuda  and  Rose 
[Ref.  llj  in  1985  in  experiments  with  R-113  condensing  on  low 
integral-fin  tubes,  confirmed  that  the  overall  heat  transfer 
coefficient  increases,  in  general,  with  decreasing  fin 
spacing.  Their  results  showed  a  600  percent  increase  in 
enhancement  over  smooth  tube  performance.  More  recent  results 
by  Marto  et  al.  [Ref.  12]  with  R-113,  demonstrated  a  700 
percent  enhancement  in  performance  over  a  smooth  tube  and  gave 
an  optimum  fin  spacing  of  0.5  mm  for  that  fluid.  Work 
performed  and  published  in  the  same  time  frame  by  Sukhatme  et 
al.  [Ref.  13]  with  R-ll  on  conventional  integral-fin  tubes  and 
special  pyramid-shaped  fin  tubes,  reported  enhancement  ratios 
of  5  to  7  for  the  low  integral  fin  tubes  and  10.3  to  12.3  for 
the  pyramid-shaped  fin  tubes. 

C.  TUBE  BUNDLE  INVESTIGATIONS 
1 .  Smooth  Tube  Bundles 

In  smooth  tube  bundles,  two  conflicting  factors  play 
a  role  in  determining  bundle  performance.  First,  the 
condensate  flowing  from  the  tubes  above  a  given  tube  in  a 
bundle  tends  to  thicken  the  condensate  film  on  that  tube, 
hence  increasing  the  resistance  to  heat  transfer.  This  effect 
is  known  as  the  condensate  inundation  effect.  Secondly, 
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droplets  from  other  tubes  striking  the  film  surface  on  a  tube 
with  a  velocity  provided  by  gravity  or  vapor  flow,  can  create 
ripples  or  waves  in  the  condensate  film  imparting  a  turbulence 
within  the  condensate  film  that  produces  an  enhancement  of 
heat  transfer  performance. 

In  1949,  Jakob  [Ref.  14]  elaborating  on  Nusselt's 
model,  predicted  that  the  vapor-side  heat  transfer  coefficient 
for  a  tube  in  a  bundle,  compared  to  the  first  tube  in  the 
bundle,  was  a  function  of  that  tube's  relative  position  in  the 
bundle.  This  model  was  based  upon  the  assumption  that  a 
continuous  laminar  sheet  of  condensate  flowing  off  the  tube 
directly  above  the  tube  considered  and  striking  the  top  of 
this  tube  further  thickened  the  condensate  film  on  this  tube. 
Jakob's  model  is  given  by: 

— ^  =  N_1^4  (2.4) 
al 

In  1958,  Kern  [Ref.  15]  proposed  a  model  based  on  the 
assumption  that  discrete  droplets  or  jets  of  fluid  from  other 
tubes  caused  ripples  in  the  condensate  film  diminishing  the 
inundation  effect.  Kern's  model  is  less  conservative  than 
Nusselt's  and,  in  many  cases,  remains  the  applied  industrial 
design  standard.  Kern's  model  is  given  by: 

^  =  N-1/6  (2.5) 
al 
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Work  by  Chen  [Ref.  16],  published  in  1961,  proposed  a  model 
that  considered  the  momentum  gain  of  falling  condensate  as 
well  as  the  condensation  of  vapor  on  sub-cooled  condensate 
droplets  or  sheets.  Chen's  model  which  is  essentially 
Nusselt's  model  times  a  factor  that  incorporates  the  phase 
change  number  is  given  by: 

=  N_1/4 [1  +  2-Ph- (N-l) ]  (2.6) 

al 

where  the  Phase  Change  number  is  given  by: 


Ph 


C  AT 
P _ 


Equation  (2.6)  is  valid  provided. 


Ph- (N-l)  <  2.0 


A  model  proposed  by  Eissenberg  [Ref.  17]  in  1972,  assumes  that 
condensate  does  not  always  drain  in  a  vertical  direction  but 
may  be  diverted  sideways  due  to  local  vapor  flow  conditions. 
In  this  case,  the  condensate  strikes  subsequent  tubes  on  their 
sides  rather  than  on  their  tops,  minimizing  inundation  effects 
on  the  condensate  film  in  the  top  portion  of  the  tube. 
Eissenberg 's  model  is  given  by: 


a 


Ji  =  .60  +  .42- (N_1/4) 
a. 


(2.7) 
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2 .  Enhanced  Tube  Bundles 

Katz  and  Geist  [Ref.  18]  in  1948,  studying  six  fin 
tubes  in  a  vertical  row,  using  R-12  among  other  working 
fluids,  found  that  the  effect  of  condensate  inundation  was 
over-predicted  by  Jakob's  model  and  proposed  that  the  exponent 
in  Equation  (2.4)  be  changed  to  0.06  for  finned  tubes. 

A  theoretical  model  was  proposed  by  Honda  et  al.  [Ref. 
19]  in  work  published  in  1987,  for  finned  tube  bundles  that 
showed  good  agreement,  within  5  to  7%,  for  data  taken 
condensing  acetone  and  R-12  on  finned  tubes.  Their  model 
relies  upon  solving  a  set  of  algebraic  equations  describing 
the  vapor-to-coolant  conjugate  heat  transfer  problem.  A 
compendium  of  the  aforementioned  models,  both  single  tube  and 
bundles,  is  provided  by  Marto  [Ref.  20]. 

As  of  the  date  of  this  writing,  the  author  is  unaware 
of  any  comprehensive  model,  substantiated  by  experimental 
results,  that  accurately  predicts  heat  transfer  performance 
for  enhanced  tube  bundles  with  varying  fin  geometries  and 
varying  pitch. 
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III.  EXPERIMENTAL  APPARATUS 

A.  CONDENSER/BOILER  TUBE  BANK  TEST  PLATFORM 
1 .  Overview 

The  composite  test  platform,  with  associated  support 
systems,  is  depicted  schematically  in  Figure  3.1.  The  boiler/ 
condenser  unit  was  fabricated  from  6.35  mm  thick,  rolled 
stainless  steel  plates  designed  to  withstand  an  absolute 
pressure  of  308  kPa.  The  top  cylinder  which  serves  as  the 
condenser  for  the  apparatus  is  1.3  0  m  in  length  with  an 
external  diameter  of  0.61  m.  The  effective  condensing  length 
is  1.22  m.  The  condensing  chamber  is  capped  on  either  end 
with  circular  stainless  steel  plates  of  0.71  m  diameter  and 
connected  to  the  chamber  with  a  flange  and  stud  assembly. 
System  integrity  is  maintained  with  a  3.20  mm  thick  rubber 
gasket.  The  end  caps  support  the  nylon  block  condenser  tube 
mounts,  stainless  steel  backing  plates,  auxiliary  condenser 
coils  (coolant  entrance  side) ,  and  mixing  chambers  (coolant 
exit  side).  The  condenser  chamber  has  five  12.7  mm  thick 
Pyrex  glass  view  ports  backed  with  12.7  mm  thick  Plexiglas. 
The  view  ports  are  127  mm  in  diameter  and  are  located  axially 
along  the  main  chamber,  angled  so  as  to  provide  top  to  bottom 
views  of  the  test  tubes  at  various  locations  along  the 
effective  condensing  length.  The  condenser  chamber  and 
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Figure  3.1  System  Schematic 


ancillary  equipment  are  depicted  graphically  and  in 
photographs  in  Figures  3.2  through  3.5. 

Internal  to  the  condenser  chamber,  a  stainless  steel 
shroud  is  fitted  over  the  effective  condensing  length  of  the 
instrumented  test  tubes  (see  Figure  3.6).  The  purpose  of  the 
shroud  is  to  channel  refrigerant  vapor  along  the  inside 
circumference  of  the  shell,  collecting  vapor  at  the  top,  and 
then  forcing  vapor  through  the  vertical  axis  of  the 
instrumented  tube  bank  into  the  shroud's  well,  where 
condensate  is  collected  and  the  remaining  vapor  is  condensed 
by  an  auxiliary  condenser.  The  auxiliary  condenser  is 
composed  of  five  helically  wound  copper  tubes  of  9.53  mm 
diameter  suspended  inside  the  shroud  well  by  cantilevered  rods 
welded  to  the  entrance  end  cap.  The  stainless  steel  shroud 
was  fabricated  with  a  glass  panel  serving  as  one  side  of  its 
stem  to  permit  viewing  of  the  instrumented  tubes  through  the 
view  ports. 

The  condenser  chamber  is  attached  to  the  boiler 
chamber  by  a  rolled  stainless-steel  cylinder  nominally  280  mm 
in  diameter  and  203  mm  in  length,  located  mid-way  along  the 
condenser  chamber's  length  allowing  condensate  to  drain  by 
gravity.  A  collar  dam  was  fitted  at  the  connecting  flange, 
for  the  initial  purpose  of  providing  a  condensate  drain  point 
to  send  condensate  to  the  auxiliary  storage  unit. 
Subsequently,  a  faster,  more-efficient  method  of  transferring 
refrigerant  was  devised  that  allowed  refrigerant  vapor  to  be 
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Fiaure  3.2  Photograph  of  Apparatus  and  Support  Systeir.s 


THERMOCOUPLE  S  1  THERMOCOUPLES  2 


17 


sent  and  condensed  in  the  auxiliary  storage  unit.  A  9.53  mm 
diameter  copper  tubing  was  provided  for  the  vapor  flow  from 
the  top  of  the  condenser  chamber  to  the  storage  tank.  A 
similar  tubing  was  provided  between  the  storage  tank  and  the 
boiling  chamber  to  allow  liquid  flow  via  gravity. 

The  boiling  chamber  is  connected  to  the  condenser  unit 
via  the  interconnecting  riser  as  previously  discussed.  The 
boiling  chamber  was  also  fabricated  from  rolled  stainless- 
steel  and  is  cylindrical  with  an  outside  diameter  of  0.61  m 
and  a  length  of  0.279  m.  The  boiling  chamber  is  fitted  with 
two  Pyrex  view  ports,  strengthened  by  Plexiglas  plates,  for 
observation  during  operation.  The  boiling  unit  is  comprised 
of  three  groups  of  tubes  (see  Figure  3.3).  In  the  simulation 
tube  bundle,  there  are  five  boiling  tubes  each  nominally  rated 
at  1.5  kW  and  located  at  the  bottom  of  the  boiler  unit.  In 
the  auxiliary  tube  bundle,  there  are  four  boiling  tubes  each 
nominally  rated  at  4  kW  and  located  two  on  either  side  of  the 
instrumented  tube  bundle.  The  instrumented  tube  bundle  is 
located  in  the  lower  center  of  the  boiling  unit  and  is 
comprised  of  35  tubes,  with  ten  active  tubes  rated  at  one  kW 
each,  five  instrumented  tubes,  and  20  dummy  tubes.  The  power 
provided  to  these  tubes  is  controlled  by  three  variac 
controllers  located  on  a  console  adjacent  to  the  test  platform 
and  fed  by  a  208  volt  power  supply.  The  variac  controllers 
are  graduated  in  one  percent  increments  of  maximum  power  for 
each  tube  bundle.  The  exact  location  of  each  tube  bundle  is 
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apparent  in  Figures  3.2  through  3.5.  It  is  a  unique  feature 
of  this  test  platform  that  both  boiling  and  condensation 
phenomena  in  a  closed  loop  system  can  be  evaluated 
simultaneously.  Details  on  the  design  and  construction  of  the 
basic  test  platform  are  available  from  Zebrowski  [Ref.  2]  and 
Murphy  [Ref.  3].  Lightoff  and  securing  procedures  for  the 
apparatus  are  listed  in  Appendix  B. 

2 .  Ancillary  Systems 

Component  equipment  that  support  the  R-114  tube  bundle 
test  platform  and  are  located  external  to  the  apparatus 
include:  (1)  an  R-114  storage  and  transfer  system,  (2)  a 
condenser  cooling  and  flow  control  system,  (3)  a  coolant 
(62.4%  by  weight  mixture  of  ethylene  glycol  and  water)  sump, 
and  (4)  an  eight-ton  refrigeration  unit. 

The  R-114  storage  and  transfer  system,  as  previously 
described,  consists  of  a  stainless-steel  cylindrical  tank 
0.350  m  in  diameter  and  0.91  m  in  length  located  on  a  rack 
above  the  coolant  sump.  Transfer  of  the  refrigerant  is 
accomplished  by  boiling  in  the  main  boiling  chamber  with  vapor 
being  sent  to  the  storage  tank  via  9.53  mm  diameter  copper 
tubing  located  in  the  top  center  of  the  condenser  chamber. 
The  vapor  is  condensed  in  the  storage  cylinder  by  means  of  a 
helical  copper  coil,  suspended  the  length  of  the  cylinder  on 
a  cantilevered  bar,  that  is  kept  cooled  by  the  water-ethylene 
glycol  mixture.  Liquid  refrigerant  can  be  returned  by  gravity 
to  the  test  platform  from  the  bottom  of  the  storage  tank 
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through  a  9.53  mm  diameter  copper  tubing  to  the  boiling 
chamber.  Faster  transfer  of  R-114  liquid  was  also  possible  if 
the  system  pressure  was  at  or  below  the  atmospheric  pressure. 
Notice  that  the  storage  tank  experiences  an  absolute  pressure 
of  about  210  kPa. 

Coolant  and  flow  control  to  the  test  apparatus  and 
ancillary  equipment  is  accomplished  by  two  different  flow  path 
systems.  Both  flow  systems  are  driven  by  two  0.5  HP  constant- 
speed  pumps  that  take  suction  on  the  main  coolant  sump.  The 
coolant  for  the  test  condenser  tube  bank  passes  from  the  pump 
discharge  through  76  mm  diameter  PVC  piping  to  a  Plexiglas 
header.  At  the  header,  flow  is  split  and  proceeds  through 
15.9  mm  diameter  Tygon  flexible  tubing  to  a  bank  of 
rotameters.  Flow  can  be  controlled  by  throttling  a  gate  valve 
located  at  the  entrance  to  each  flow  meter.  Coolant  flow 
leaves  the  exit  of  each  flow  meter  and  proceeds  to  the 
instrumented  condenser  tube  bank  through  flexible  Tygon 
tubing.  At  the  exit  of  the  main  condenser  chamber,  coolant 
leaves  each  tube  and  flows  through  flexible  Tygon  tubing  to 
individual  mixing  chambers  (Figure  3.7).  From  the  mixing 
chambers,  coolant  passes  through  flexible  Tygon  tubing  to  a 
central  Plexiglas  header  suspended  above  the  test  platform, 
and  the  collective  coolant  is  piped  back  to  the  main  sump 
through  76  mm  diameter  PVC  pipe.  Auxiliary  coolant  flow 
leaves  the  pump  discharge  through  PVC  piping  and  proceeds  to 
a  two-way  ball  valve,  where  the  flow  is  either  sent  to  the 
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Figure  3.7  Mixing  Chamber  Schematic 
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R-114  storage  tank  and  back  to  the  main  sump  or  is  sent 
through  PVC  pipe  to  a  large  rotameter.  From  the  auxiliary 
flowmeter,  coolant  enters  a  central  header  and  is  then  split 
to  flow  to  the  five  auxiliary  condenser  coils  and  back  to  the 
coolant  sump.  Flow  through  each  auxiliary  condenser  coil  is 
controlled  by  a  valve  located  at  the  coil  inlet.  The  coolant 
sump  has  a  1.81  cubic  meter  capacity  and  is  constructed  of 
12.7  mm  thick  sheets  of  Plexiglas.  The  coolant  is 
approximately  62.4%  ethylene  glycol  and  37.6%  water  (by 
weight) . 

The  coolant  is  chilled  by  an  externally-located,  8  ton 
refrigeration  system  that  continually  re-circulates  sump 
coolant  with  a  0.75  HP  pump.  The  refrigeration  system  is 
capable  of  maintaining  a  sump  temperature  between  -21  C  and 
ambient  temperature. 

3 •  Instrumentation 

The  coolant  temperature  r.  se  is  measured  by  series- 
connected  thermopiles  having  ten  junctions  on  either  end. 
These  thermopiles  were  fabricated  using  type-T  thermocouple 
wire.  The  ends  of  the  thermopiles  are  inserted  in  stainless 
steel  wells  with  copper-plugged  tips  located  at  the  entrance 
to  each  individual  condenser  tube  and  at  the  exit  of  this 
tube's  mixing  chamber.  Inlet  coolant  temperature  is  measured 
with  type-T  thermocouples  also  located  in  the  entrance  wells. 
Refrigerant  liquid  and  vapor  temperatures  are  measured  with 
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type-T  thermocouples  inserted  in  stainless-steel  wells  at  the 
locations  indicated  in  Figures  3.2  through  3.5. 


System  pressure  is  monitored  through  a  calibrated 
pressure-vacuum  gage  valid  over  a  range  of  30  inches  mercury 
(vacuum)  to  30  pounds  per  square  inch  (gage  pressure) . 

4 .  Tube  Bundle  Data  Acquisition  and  Reduction 

A  Hewlett-Packard  9S16A  computer  was  used  to  control 
a  Hewlett-Packard  3497A  Automatic  Data  Acquisition  System, 
which  read  the  output  of  the  thermopiles  and  thermocouples. 
Readings  were  made  in  millivolts  and  were  converted  to 
temperature  readings  in  the  data  reduction  program.  The 
channels  read  by  the  data  acquisition  system  are  listed  in 
Table  3.1. 

5 .  Tubes  Tested 

Two  sets  of  tubes  were  tested.  The  first  was  a  smooth 
copper  tube  (inside  diameter  13.26  mm,  outside  diameter  15.88 
mm)  and  the  second  was  a  low  integral-fin  copper-nickel  tube 
(inside  diameter  10.16  mm,  root  diameter  14.00  mm,  outside 
diameter  15.88).  These  two  sets  of  tubes  were  tested,  in 
bundles  and  individually,  with  R-114  and  R-113. 
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TABLE  3.1 


CHANNEL  ASSIGNMENTS  ON  DATA  ACQUISITION  SYSTEM 


Channel  Numbers 


Measurement 


20 

21 

22 

23 


Thermopile  1st  Tube 
Thermopile  2nd  Tube 
Thermopile  3rd  Tube 
Thermopile  4th  Tube 
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IV.  SYSTEM  OPERATION  AND  DATA  REDUCTION 

A.  TEST  PLAN  AND  MODIFICATION  OF  DESIGN  FOR  CONDENSER 

TEST  APPARATUS 

1 .  Ability  to  Change  Working  Fluids 

Although  envisioned  in  the  original  apparatus  design 
by  Zebrowski  [Ref.  2]  and  Murphy  [Ref.  3],  no  physical  system 
existed  for  the  change-out  and  storage  of  working  fluids.  The 
system  was  designed  as  a  general  test  platform  to  evaluate  the 
performance  of  various  refrigerants  and  steam  during 
condensation  on  a  variety  of  test  tubes. 

The  R-114  storage  and  transfer  system  as  described  in 
Chapter  III,  ancillary  systems,  was  designed  and  built  to 
allow  for  the  storage  of  refrigerants  that  evaporate  at 
atmospheric  temperature  and  pressure.  This  additional 
capability  not  only  conserves  expensive  refrigerants  but 
facilitates  speedy  change  of  these  working  fluids  during  the 
evaluation  of  a  particular  tube  with  a  variety  of 
refrigerants. 

2 .  Tube  Alignment  and  Reduction  of  Tubes  in  Test  Bundle 

Visual  inspection,  during  preliminary  experimental 

runs,  revealed  that  bowing  of  the  tube  bundle  was 
detrimentally  affecting  condensation  patterns.  As  a  result  of 
misalignment  of  the  endcaps  during  fabrication,  bending 
moments  of  varying  magnitudes  and  directions  resulted  in 
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condensate  flow  striking  subsequent  tubes  in  the  vertically- 
oriented  bundle  at  different  positions  around  the 
circumference  of  the  tubes.  Upon  disassembly  of  the  end 
plates  of  the  apparatus  and  comparison  with  the  proposed 
drawings  for  the  DDG-51  refrigeration  condenser  provided  by 
the  David  Taylor  Research  Center,  miscalculations  in  the  tube 
pitch  were  revealed.  The  vertical  pitch  required  for  the 
project  is  35.74  mm,  centerline  to  centerline.  The  nylon 
bundle  plates,  as  manufactured,  had  a  pitch  that  was 
considerably  less.  In  an  effort  to  correct  the  tube  mis¬ 
alignment  problem,  a  low-powered  laser  was  proposed  to  align 
the  bundle  tubes.  A  new  nylon  endplate  was  fabricated  for  the 
exit  endcap,  at  the  correct  pitch  with  the  tube  bundle  reduced 
to  four  vertical  tubes, due  to  constraints  in  the  endcap 
openings.  An  aluminum  mount  for  the  laser  was  fabricated  by 
machine  shop  personnel  that  fit  snugly  into  the  tube 
penetrations  of  the  nylon  block.  Several  test  projections 
were  conducted  on  a  plastic  template  fitted  to  the  entrance 
endcap  to  minimize  misalignment  possibilities.  The  plastic 
template  was  scribed  and  cut  with  particular  attention  paid  to 
correct  pitch  and  used  to  manufacture  the  entrance  nylon 
endplate.  Subsequent  visual  inspection  during  experimental 
runs  revealed  that  misalignments  and  bowing  had  been 
corrected. 
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3 .  Vapor  Superheat  Problems 


During  initial  experimental  runs,  observed 
discrepancies  in  the  measured  temperature  increases  across  the 
condenser  tube  lengths  when  compared  to  predicted  temperature 
increases,  prompted  consideration  of  the  possibility  of  vapor 
superheat  occurring  as  vapor  (at  saturation  temperature  less 
than  room  temperature)  flowed  from  the  boiling  chamber  through 
the  riser  and  up  the  sides  of  the  condensing  chamber.  A  two¬ 
pronged  strategy  was  developed  to  minimize  this  possible 
effect.  First,  the  apparatus  shell  was  completely  insulated 
with  12.7  mm  thick  foam  insulation  and  all  ancillary  tygon 
tubing  was  insulated  with  double-wrapped,  3.18  mm  thick  foam 
insulation.  Secondly,  all  experimental  runs  were  performed  at 
vapor  saturation  temperatures  above  the  ambient  temperature. 

4 .  Contamination  Problems 

Visual  inspection  of  the  condensing  tubes  during 
operation,  correlated  with  observed  disparities  in  measured 
temperature  increases  across  condenser  tubes,  revealed  a 
structured,  crystalline  surface  or  matrix  formation  on  the 
tube  surfaces  that  inhibited  heat  transfer.  The  first 
appearance  of  the  inhibiting  matrix  occurred  during  runs  with 
smooth  copper  tubes  and  with  the  coolant  inlet  temperature  at 
approximately  -20  C.  The  working  fluid,  R-114,  was 
transferred  to  the  storage  unit  except  for  an  oily  residue 
that  remained  in  the  bottom  of  the  boiling  chamber.  This 
residue  was  not  analyzed,  but  was  assumed  to  be  machining  oil 
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that  remained  after  apparatus  fabrication.  The  system  was 
flushed  twice,  once  with  acetone  and  once  with  R-113,  and  the 
boiling  chamber  was  scrubbed  out.  The  smooth  copper  tubes 
were  re-installed  after  a  complete  cleaning  of  their  outside 
surfaces,  and  the  experimental  runs  were  repeated  under  the 
same  conditions.  Visual  inspection  again  revealed  the 


presence  of 

the 

matrix 

formation  on 

the 

tubes,  and  the 

possibility 

of 

water 

contamination 

in 

the  R-114  was 

conjectured  as  a  probable  cause.  The  matrix  was  thought  to  be 
ice  crystals  solidified  on  the  tube  surface.  The  decision  was 
made  to  continue  experimental  runs  with  a  different 
refrigerant  (R-113)  and  a  different  type  of  condensing  tube 
(copper-nickel  1024  fpm  low  integral-fin  tube) ,  with  the  aim 
of  understanding  the  conditions  that  caused  the  phenomena  to 
exist.  First,  a  run  was  made  using  R-114  with  the  copper- 
nickel  finned  tube  under  exactly  the  same  conditions  as 
described  in  the  smooth  copper  tube  runs.  Visual  inspection, 
during  this  run,  revealed  an  apparent  thickening  of  the  fins 
at  the  top  of  the  tube  when  compared  visually  with  the  second 
tube  (see  Figure  4.1(a)  fin  normal  appearance  and  Figure 
4.1(b)  fin's  thickened  appearance).  This  appearance  of 
thickening  gave  credence  to  the  belief  that  water 
contamination  was  causing  an  ice  layer  to  form  on  the  top  tube 
between  fins.  The  refrigerant  was  then  changed  to  R-113  and 
the  run  was  repeated,  with  the  inlet  cooling  temperature  at 
approximately  -20  C.  Again,  a  marked  thickening  of  the  finned 
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surface  appeared  on  the  top  of  the  first  tube.  The  next 
experimental  run  was  conducted  on  the  same  finned  tube  with  R- 
113,  but  at  a  coolant  inlet  temperature  slightly  above  0  C. 
In  this  situation,  no  visible  thickening  of  the  finned  surface 
was  observed,  but  disparities  in  temperature  increases  across 
tube  lengths  persisted.  The  run  was  repeated  several  days 
later  to  observe  whether  the  phenomenon  consistently  repeated 
itself.  Finally,  the  smooth  copper  tube  was  substituted  for 
the  finned  tube  and  the  experimental  run  repeated  at  the 
higher  coolant  inlet  temperature.  Noticeable  surface 
contamination,  having  an  orange  peel  appearance,  recurred.  In 
addition,  during  this  run,  there  were  indications  that  a 
pocket  of  non-condensable  gases  formed  in  the  top  of  the 
condensing  chamber.  As  every  effort  had  been  made  to  evacuate 
the  apparatus  prior  to  commencing  the  experimental  run,  and  no 
evidence  exists  to  support  outside  leakage  of  air  into  the 
system,  the  possibility  exists  that  contamination  within  the 
system  produced  the  non-condensable  gases  during  boiling. 

Two  possible  sources  of  the  contamination  are 
conjectured.  The  first  is  that  high  temperatures  on  the 
heater  tubes  and  poor  circulation  in  the  boiling  chamber  are 
combining  to  break  down  the  refrigerant  molecules  producing 
some  type  of  hydro-carbon.  The  second  is  that  a  chemical 
reaction  is  occurring  between  the  refrigerants  and  gasket 
material  in  the  apparatus  producing  a  hydro-carbon.  Time  and 
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funding  limitations  have  prevented  further  investigation  of 
the  problem. 

5 .  Summary 

The  original  test  plan  called  for  the  testing  of 
various  enhanced  surfaces  in  a  simulated  bundle  during 
condensation  of  R-114.  The  contamination  problem  described 
above,  coupled  with  the  time  delays  inherent  in  achieving 
solutions  to  the  other  aforementioned  encountered  problems, 
severely  limited  the  results  of  this  thesis. 

B.  DATA  REDUCTION 

1 •  Description  of  Program  Capabilities 

The  computer  program,  DRP1F,  that  collects  and 
processes  raw  data,  in  conjunction  with  a  Hewlett-Packard 
computer/data  acquisition  system  described  in  Chapter  II,  is 
listed  in  its  entirety  in  Appendix  C.  The  program  is  designed 
to  calculate  and  plot  heat  transfer  parameters  for  a  variety 
of  different  tube  bundles,  utilizing  R-114,  R-113,  or  steam  as 
the  working  fluid.  The  program  has  the  added  capability  of 
allowing  testing  for  single  tube  performance.  The  program 
consists  of  five  main  sections,  as  follows: 

1.  Driver  Program, 

2.  Main  Program, 

3.  Property  Subroutines, 

4.  Modified  Wilson  Plot  Subroutine. 

5.  Plotting  Subroutines. 
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Of  these  five  sections,  only  the  modified  Wilson  Plot 
subroutine  will  be  described  in  lengthy  detail. 

The  driver  program  permits  the  user  to  take  data, 
reprocess  data,  or  plot  re-processed  data  through  various 
subroutines.  ine  driver  program  is  listed  in  lines  1000 
through  1125,  in  Appendix  C. 

The  main  subroutine  (lines  1130-2315)  can  be  divided 
into  five  parts.  The  first  part  allows  the  user  to  select  the 
physical  parameters  used  in  data  reduction.  The  selection  of 
parameters  consists  of  the  working  fluid  to  be  used  (R-114,  R- 
113,  or  steam),  the  vapor  saturation  temperature  to  be  used 
(derived  from  averaged  thermocouple  readings  in  either  the 
vapor  section  or  liquid  region  of  the  apparatus) ,  the 
instrumented  test  tube  type,  and  whether  data  will  be  taken 
on  a  bundle  or  individual  tubes.  The  second  part  allows  the 
user  to  reprocess  data  with  the  same  selection  of  physical 
parameters  described  above,  but  calls  the  modified  Wilson  Plot 
subroutine  to  calculate  the  inside  and  outside  coefficients 
used  in  the  correlations.  Basic  data  reduction  takes  place  in 
the  third  part  (lines  1945-2125) .  Subroutines  are  called  to 
calculate  the  properties  of  the  ethylene  glycol-water 
solution.  These  properties,  with  the  physical  tube 
dimensions,  are  used  to  calculate  velocities  and  Reynolds 
numbers.  Low  coolant  velocities  and  high  viscosities  prompted 
the  use  of  twisted  tape  inserts  (thickness  0.559  mm,  with  a 
pitch  for  a  180  degree  twist  of  three  times  the  tube's  inner 
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diameter)  to  increase  the  inside  heat  flux.  The  inside 
Nusselt  number  with  a  twisted  tape  insert,  provided  by  Hong 
and  Bergles  [Ref.  21],  is  given  by  the  correlation: 


Re 


H5.  =  .  5.17211  +5.4838-10-3-  (Pr0'7)  •  ^ V* 


(4.1) 


where  the  Reynolds  number  for  coolant  flow  was  given  by: 


Re 

s 


4*r/ir*ri*  (Dj_  -  46) 


(4.2) 


where,  6  is  tape  thickness. 

The  outside  heat  transfer  coefficient  is  then  given  by 
the  well-known  summation  of  resistances  to  heat  transfer: 


UX-“=t-+Rf  +  Rm  +  ^ 


(4.3) 


o  o  a  A 
o  o 
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which  algebraically  reduces  to: 
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a  = 
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rf-  - R  A  -R-A 
U  mo  r  o  —  A. 


(4.4) 
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The  heat  transferred  to  the  coolant  (Q)  is  given  by  the 
relationship: 


Q  = 


r-C  • (AT) 

p 


(4.5) 
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The  heat  flux  (QM )  is  subsequently  calculated  by  dividing  by 
the  outside  surface  area: 


Q" 


(4.6) 


where  A0  =  ^do  -L.  And,  the  overall  heat  transfer  coefficient 
(U0)  is  given  by: 


U 


o 


=  Q" 

LMTD 


(4.7) 


where  the  Log  Mean  Temperature  Difference  (LMTD)  is  defined 
to  be: 


IMID  = 


AT 


T  -Tc- 
,  r  sat  1, 

log  [m  _l rn-,  J 

sat  o 


(4.8) 


It  should  be  noted  that  the  wall  resistance  due  to  fouling  was 
assumed  to  be  negligible  for  the  purposes  of  calculation. 

The  fourth  part  of  the  main  subroutine  creates  a  raw 
data  file  (lines  1680-1835) ,  allowing  subsequent  reprocessing 
by  the  modified  Wilson  subroutine.  The  fifth  part  of  the  main 
subroutine  provides  a  printed  output  both  while  taking  initial 
data  and  subsequently  after  reprocessing  data. 

The  third  major  section  of  the  computer  program 
calculates  fluid  properties  through  called  functions  for  both 
the  working  fluid  and  the  coolant.  The  calculated  coolant 
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properties  as  a  function  of  temperature  (lines  2330-2615) ,  are 
kinematic  viscosity,  specific  heat,  density,  Prandtl  number, 
and  conductivity. 

The  fifth  major  section  of  the  program  provides 
plotting  routines  for  the  output  files  generated  in  the 
program.  The  relationships  graphically  displayed  are  heat- 
transfer  coefficient  «ratios  (either  based  on  the  first  tube  in 
the  bundle  or  as  a  ratio  of  the  Nusselt  value)  plotted  against 
tube  position,  heat  transfer  coefficient  plotted  against 
either  the  heat  flux  or  temperature  rise  of  the  coolant,  and 
the  X-Y  plot  generated  from  the  modified  Wilson  Plot  results. 

2 .  Modified  Wilson  Plot 

The  modified  Wilson  Plot,  as  outlined  by  Marto  [Ref. 
22],  accomplishes  an  indirect  measurement  of  the  outside  heat- 
transfer  coefficient.  The  implicit  assumption  is  that  the 
overall  heat  transfer  coefficient  (UG)  is  reliably  known  from 
data  and,  therefore  a  summation  of  heat  transfer  resistances 
is  assumed.  This  summation  relationship  is  given  by: 


r  -  =■ +  Vo +  Vo +  =r<5r> 


(4.9) 


“i  1 


where  resistance  due  to  wall  fouling  is  assumed  equal  to  zero. 
The  summation  equation  is  transformed  to  a  linear 
relationship,  as  follows: 


kr-Vo]  = 


Jl  +  ± 

VA.  - 

l  a.  a 

1  o 


(4.10) 
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where : 


a. 

1 


=  C.  (— ) 
l 'd.  ' 
1 


a 


=  CoF 


(4.11) 


which  results  in  the  simple  linear  form,  Y  =  mX  +  b.  It 
should  be  noted  that  Theta  (defined  in  line  3290)  is  derived 
from  the  Hong  and  Bergles  [Ref.  20]  relationship  for  the 
inside  Nusselt  number,  given  by: 

Nu.  =  =  C.  [1  +5.4838-10_3(Pr°*7)  •  (— : -) 1  “ 25] 1/2  (4.12) 

i  a  x  y 

hence. 
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[1  +  5.4838*10-3(Pr°'7)  •  (_®)  1-25]1/2 


(4.13) 


Further,  F  (defined  in  line  3330),  is  derived  from  Nusselt’s 
relationship  for  a  horizontal  cylinder  subjected  to  a  constant 
heat  flux  [Ref.  4],  given  by: 


=  .655  [A3,  p2*g*Ah yn  ’  Dq  •  Q"  ] 1/3 


(4.14) 


Hence, 


F  =  [A3-p2-g*Ahv/n-Do-Q"]1/3 


(4.15) 
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X  and  Y  values  are  calculated  from  raw  data  and  the  data  are 
fit  with  a  least  squares  approximation.  Initial  assumed 
values  are  taken  from  the  aforementioned  correlations,  and  the 
solutions  iterated  to  find  the  inside  coefficient  (CL)  and  the 
outside  coefficient  (C0)  that  fit  the  data,  where: 


D  -F 

Y  =  [rf  -  R  (A  )]F  ;  X  =  -~- 
U  m  o  A 

o 


the  slope  (m)  is  given  by: 


(4.16) 


m  - 


C. 

1 


and  the  intercept  (b)  is  given  by: 


The  accuracy  of  this  method  relies  heavily  upon  the  number  of 
data  points  taken  and  the  range  of  velocities  utilized. 
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V.  RESULTS  AND  DISCUSSION 

A.  EXPERIMENTAL  SEQUENCE 

The  data  runs  are  summarized  in  Table  5.1  according  to 
tube  type,  refrigerant  used,  and  approximate  operating 
parameters.  Specifics  for  each  run  are  provided  in  the 
following  section  under  the  run  title. 


TABLE  5.1 

SEQUENTIAL  LISTING  OF  DATA  RUNS 


Run 

Title 

Tube 

Tvoe 

Refrigerant 

Used 

Coolant 
T  inlet 

Vapor 

T  saturation 

SMT02 

smooth 

R-114 

-20  C 

17.8 

C 

CNFT01 

finned 

R-114 

-21  C 

38.3 

C 

CNFT02 

finned 

R-113 

-21  C 

64.4 

C 

CNFT03 

finned 

R-113 

-3  to  +6 

C 

52.2 

C 

CNFT04 

finned 

R-113 

-3  to  +6 

c 

46.3 

C 

SMT03 

smooth 

R-113 

-1  to  +6 

c 

56.7 

C 

B.  EXPERIMENTAL  RESULTS 

1 .  Smooth  Tubes  with  R-114  and  Low  Coolant  Temperature 
SMT02  was  an  experimental  run  made  with  four  smooth 
copper  tubes  while  condensing  R-114.  Vapor  saturation 
temperature  was  17.8  C.  Ambient  temperature  was  20  C.  The 
ethylene  glycol  coolant  inlet  temperature  was  -20  C.  The 
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experiment  was  conducted  following  a  cleaning  of  the  boiling 
chamber,  two  flushings  of  the  apparatus,  once  with  acetone  and 
once  with  R-113,  evacuation  to  27.5  in  Hg.,  and  filling  the 
system  from  an  unopened  cylinder  of  R-114.  After  filling  the 
system,  the  apparatus  pressure  stabilized  at  17  psig.  The 
system  pressure  was  lowered  to  10  psig  by  opening  the  coolant 
flow  through  the  auxiliary  condenser.  When  the  desired  gage 
pressure  was  reached,  the  lowest  flow  rate  was  set  through  the 
instrumented  tube  bundle,  and  all  three  heater  units  were  set 
for  power  levels  corresponding  to  approximately  10  kW  total. 
Gage  pressure  was  maintained  nearly  constant  by  controlling 
coolant  flow  through  the  auxiliary  condenser.  The  system 
required  approximately  35  minutes  from  lightoff  to  reach  a 
steady  state  condition,  indicated  by  nearly  constant  gage 
pressure.  The  R-114  appeared  clear  at  the  commencement  of  the 
data  run.  By  the  time  the  apparatus  reached  a  steady  state 
condition,  the  contamination  matrix  was  fully  formed  and 
visible  on  the  top  two  tubes  in  the  bundle.  Data  were  taken 
at  nine  different  coolant  velocities,  after  the  voltage 
indicator  on  the  data  acquisition  system  for  the  bottom  tube 
thermopile  appeared  to  reach  a  fixed  value  following  each 
velocity  change.  Data  for  the  tube  bundle  run  (Figure  5.1) 
clearly  demonstrates  an  observed  uncertainty  of  about  7%  in 
the  overall  heat  transfer  coefficient  (U0)  at  velocities  close 
to  0.9  m/s.  This  observed  uncertainty,  in  contrast  to  the 
maximum  calculated  uncertainty  of  4.5%  (see  Appendix  A), 
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aure  5.1  SitT02 — Tube  Bundle  Performance 


indicates  the  effect  of  the  contamination  on  the  measurement 
of  the  overall  heat  transfer.  The  most  significant  result  of 
the  contamination  is  the  addition  of  an  un-quantif iable 
resistance  to  heat-transfer  that  has  produced  heat-transfer 
coefficients  approximately  half  of  those  predicted  by  the 
Nusselt  correlation.  Further,  it  is  apparent  that  during 
bundle  operation,  the  contamination  degrades  the  performance 
of  the  tubes  in  a  graduated  manner  with  the  first  tube  being 
the  most  adversely  affected  and  the  third  tube  appearing  the 
least  affected. 

Upon  completion  of  the  bundle  run,  the  flow  in  all  of 
the  instrumented  tubes  was  shut  off.  Gage  pressure  was 
maintained  at  10  psig  by  controlling  coolant  flow  through  the 
auxiliary  condenser.  The  coolant  inlet  temperature  remained 
approximately  -20  C.  The  apparatus  remained  in  this  condition 
until  the  tube  surfaces  appeared  to  dry  off,  which  took 
normally  about  ten  minutes.  Data  were  then  taken  at  the  same 
nine  flow  rates,  as  for  the  bundle,  but  with  flow  through  only 
one  instrumented  tube  at  a  time,  starting  with  the  top  tube. 
Between  velocity  changes,  approximately  10  minutes  was  allowed 
for  the  temperature  changes  in  the  tube  to  take  effect.  Data 
runs  were  made  only  after  inundation  from  the  previously 
tested  tube  was  no  longer  present.  The  R-114  in  the  boiling 
chamber  remained  visually  clear  for  all  single  tube  runs. 
Again,  as  in  the  bundle  data,  data  from  the  single  tube  runs 
(Figure  5.2)  fall  within  an  acceptable  uncertainty  band. 
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From  the  conditions  under  which  these  runs  were  made,  it  is 
noteworthy  and  logical  that  the  effect  of  the  contamination 
increases  in  magnitude  as  a  run  progresses.  The  first  tube 
gives  the  lowest  results  possibly  because  of  residual 
contamination  from  the  bundle  run.  The  first  and  second  tubes 
produced  approximately  the  same  results  as  previously 
discussed  in  the  bundle  run,  however  the  third  and  fourth 
tubes  gave  lower  results  when  operating  as  single  tubes  than 
when  operating  in  a  bundle,  possibly  because  condensate 
inundation  that  provides  a  rinsing  effect  in  the  bundle 
operation, is  no  longer  present  during  single  tube  operation. 

2 .  Finned  Tubes  with  R-114  and  Low  Coolant  Temperature 
CNFT01  was  an  experimental  run  made  with  copper-nickel 
low  integral-fin  tubes  condensing  R-114.  Vapor  saturation 
temperature  was  18.3  C.  The  ethylene  glycol  coolant  inlet 
temperature  was  at  -21  C.  The  experiment  commenced  following 
the  completion  of  the  SMT02  single  tube  runs,  transfer  of  the 
R-114  to  the  auxiliary  storage  unit,  installation  of  the 
copper-nickel  finned  tubes,  evacuation  of  the  apparatus  to 
27.5  in  Hg,  and  a  refill  of  R-114  from  the  storage  unit.  This 
process  took  approximately  six  hours  to  complete.  The  same 
data  taking  procedure  as  outlined  for  the  smooth  tube  runs  waL, 
followed.  The  appearance  of  the  R-114  was  clear  at  the 
commencement  of  the  run.  By  the  time  the  system  reached 

steady  state,  a  marked  thickening  of  the  fin  areas  on  the  top 
tube  when  compared  to  the  fins  of  the  second  tube  was  apparent 


(see  Figure  4.1).  The  R-114  appeared  clear  at  the  end  of  the 
bundle  run.  Data  from  the  bundle  run  (Figure  5.3)  show  an 
increase  in  the  overall  heat  transfer  coefficient  of 
approximately  100%  when  compared  to  the  bundle  results  for  the 
smooth  tubes  SMT02  (Figure  5.1).  The  coolant  velocity 
differences  between  the  smooth  tubes  and  the  finned  tubes  is 
due  to  the  smaller  inside  diameter  of  the  finned  tubes  at  the 
prescribed  mass  flow  rate.  There  is  no  disparity  between 
calculated  and  observed  uncertainties  given  the  calculated 
uncertainty  band,  but  performance  remains  approximately  half 
of  that  expected  when  compared  to  the  values  obtained  with  the 
Nusselt  correlation  and  enhancement  ratios  reported  by  other 
investigations.  Once  again, the  first  and  second  tubes 
demonstrate  the  greatest  degradation  in  performance  due  to  the 
contamination.  The  third  and  fourth  tubes  either  receive  less 
contamination  or  derive  enhancement  from  condensate  inundation 
flow.  From  observations,  the  unknown  contaminant  had  a  higher 
surface  tension  than  R-114,  but  to  what  extent  the  contaminant 
kept  fin  root  areas  flooded  and  subsequently  negated  the 
enhanced  surface  effect,  remains  obscured  in  the  uncertainties 
of  the  measurements. 

The  single  tube  runs  followed  the  procedures  outlined 
for  the  single  tube  runs  of  the  smooth  copper  tubes.  However 
at  the  end  of  each  single  tube  run,  a  marked  thickening  of  the 
fin  surface  at  the  top  of  the  tube  was  noticeable.  Vapor 
saturation  temperature  and  coolant  inlet  temperature  were  the 
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0.95  1.10  1.25  1.10  1.55  1.70  1.05 


same  as  in  the  bundle  run.  Data  from  the  single  tube  runs 
(Figure  5.4)  demonstrate  an  enhancement  of  approximately  76% 
in  the  overall  heat  transfer  coefficient  when  compared  to  its 
counterpart  in  SMT02 .  The  spread  in  the  data  is  approximately 
7%  and,  given  the  calculated  uncertainty  band,  appears 
reasonable.  Performance  remains  lower  than  expected  from 
existing  studies,  and  clearly  the  magnitude  of  the 
contamination  effect  increases  as  each  run  proceeded. 

3 .  Finned  Tubes  with  R-113  and  Low  Coolant  Temperature 
Upon  completion  of  the  single  tube  runs  described  in 
CNFT01,  R-114  was  transferred  to  the  storage  unit.  The 
apparatus  was  drained  of  residual  refrigerant,  the  system  was 
evacuated  to  27.5  in  Hg,  and  filled  with  freshly  distilled  R- 
113  by  using  apparatus  vacuum  to  promote  flow  from  the  R-113 
container.  After  filling,  system  vacuum  was  at  11  inches  Hg. 
The  system  was  brought  to  a  gage  pressure  of  5  psig 
corresponding  to  a  vapor  saturation  temperature  of  64.4  C. 
Coolant  inlet  temperature  remained  the  same  as  described  in 
CNFTOl,  -21  C.  Vapor  saturation  temperature  was  monitored 
through  the  system  pressure,  and  controlled  by  regulating 
coolant  flow  through  the  auxiliary  condenser.  The  same  nine 
data  points  described  previously  were  taken  in  accordance  with 
the  stated  procedure.  The  appearance  of  the  R-113  at  the 
beginning  of  the  run  was  clear  and  clean.  Upon  the  conclusion 
of  the  run,  however,  the  R-113  had  a  slight  yellowish  tinge. 
During  the  run,  the  same  marked  thickening  of  the  fin  surface 


48 


described  in  CNFT01,  and  shown  schematically  in  Figure  4.1 
occurred  on  the  top  tube.  The  data  from  the  bundle  run 
(Figure  5.5)  demonstrates  only  a  slight  improvement  in  the 
overall  heat  transfer  coefficient  when  compared  to  the  bundle 
data  for  SMT02  (Figure  5.1),  and  in  fact  represents  an 
approximate  80%  decrease  in  the  overall  heat  transfer 
coefficient  when  compared  to  the  performance  of  the  bundle  in 
CNFT01  (Figure  5.3).  In  view  of  the  fact  that  these  data  were 
collected  after  transferring  the  R-114  to  the  storage  unit  and 
filling  with  freshly  distilled  R-113,  it  can  only  be  assumed 
that  the  contamination  was  freshly  created  with  the  R-113  and, 
in  fact  is  associated  with  the  apparatus.  It  is  also 
possible,  based  on  these  results,  to  infer  the  presence  of 
non-condensable  gases  generated  by  the  contamination,  although 
no  indications  of  non-condensable  gases  were  detected  during 
the  run.  In  view  of  these  facts,  the  reliability  of  the 
calculated  values  for  the  overall  heat  transfer  coefficient  is 
extremely  doubtful  even  though  data  spread  falls  within  the 
acceptable  uncertainty  band. 

Upon  completion  of  the  bundle  run,  data  were  taken  for 
each  of  the  tubes  individually  at  the  same  system  conditions 
described  for  the  bundle  run,  and  in  accordance  with  the 
procedure  described  for  single  tube  runs.  The  appearance  of 
the  R-113  at  the  completion  of  the  single  tube  runs  was 
slightly-darkened,  with  a  yellow  tinge,  but  otherwise  appeared 
free  from  any  floating  contaminants.  On  each  of  the  single 
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tube  runs,  a  marked  thickening  of  the  fin  surface  developed  in 
the  course  of  the  run.  There  exists  reasonably  good  agreement 
between  the  single  tube  data  (Figure  5.6)  and  the  single  tube 
data  from  CNFT01  in  the  calculation  of  the  overall  heat 
transfer  coefficient,  and  no  disparity  exists  between  the  data 
spread  and  the  calculated  uncertainty  band.  The  apparatus  was 
evacuated  for  five  minutes  prior  to  beginning  the  single  tube 
runs  and  if  non-condensable  gases  were  present  during  the 
bundle  run,  this  explains  the  improved  performance  of  the 
tubes  during  individual  runs. 

4 .  Finned  Tubes  with  R-113  and  High  Coolant  Temperature 
This  data  run  was  made  with  the  same  copper-nickel 
tubes  described  in  CNFT01  and  CNFT02.  The  system  was  allowed 
to  sit  overnight  at  11  in  Hg,  while  the  sump  temperature  was 
brought  to  an  operating  range  of  -3  C  to  +6  C.  No  change  in 
the  system  vacuum  over  this  12  hour  period  occurred.  The 
decision  to  run  the  system  at  a  higher  coolant  temperature  was 
motivated  by  a  belief  that  the  contamination  of  the  system  was 
water,  and  that  at  the  colder  coolant  temperature,  an  ice 
sheath  was  forming  on  the  tubes  inhibiting  heat  transfer.  The 
appearance  of  the  R-113  before  commencement  of  the  bundle  run 
was  a  medium  yellow  tinge.  The  system  was  brought  to  a  steady 
state  pressure  corresponding  to  a  vapor  saturation  temperature 
of  52.2  C.  Difficulty  was  encountered  in  maintaining  positive 
system  pressure  with  the  auxiliary  condenser  operating.  The 
auxiliary  condenser  was  therefore  turned  off  and  steady  state 
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VEL.  (M/S) 

CNFT02  —  Individual  Tube  Performance 


was  maintained  by  slight  power  adjustments  to  the  heating 
elements  in  the  boiling  chamber.  Nine  sets  of  data  were  taken 
in  accordance  with  the  procedure  for  bundle  runs  previously 
described.  No  observable  thickening  of  fin  surfaces  during 
the  run  occurred.  The  data  from  the  bundle  run  (Figure  5.7) 
shows  a  slight  increase  in  overall  heat  transfer  coefficient 
over  the  bundle  results  obtained  in  CNFT02  (Figure  5.5)  and 
this  is  possibly  explained  by  a  decrease  in  the  viscosity  of 
the  contaminant  when  exposed  to  the  warmer  tube  surface. 

Single  tube  runs  were  made  on  each  of  the  tubes,  in 
accordance  with  procedures  outlined  previously,  at  the  same 
coolant  inlet  temperature  as  in  the  bundle  run,  but  at  a 
system  pressure  corresponding  to  a  vapor  saturation 
temperature  of  64.4  C.  No  change  in  the  appearance  of  the  R- 
113  was  detected.  Again  no  thickening  of  the  fin  surfaces  was 
detected.  The  single  tube  results  (Figure  5.8)  demonstrate  a 
corresponding  range  of  values  in  the  overall  heat  transfer 
coefficient  when  compared  to  the  results  for  single  tubes  in 
CNFT02  (Figure  5.6),  but  is  significantly  higher  than  the 
bundle  results.  It  is  possible  that  non-condensable  gases 
were  present  in  bundle  operation,  but  were  eliminated  by 
evacuation  prior  to  single  tube  data  being  taken.  No 
indications  of  non-condensable  gases  were  detected  during 
bundle  or  single  tube  operation. 
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VEL.  (M/S) 

Ficrure  5.7 — CNFT03 — Tube  Bundle  Performance 


VEL.  (M/S) 

-CNFT03 — Individual  Tube  Performance 


5 .  Finned  Tubes  with  R-113  and  High  Coolant  Temperature 

This  data  run  was  made  in  order  to  see  if  the  results 
produced  in  CNFT03  were  repeatable.  The  system  operating 
pressure  corresponded  to  a  vapor  saturation  temperature  of 
46.3  C,  and  the  coolant  inlet  temperature  remained  unchanged 
from  the  coolant  temperature  in  CNFT03.  Nine  data  sets  were 
taken  in  accordance  with  the  bundle  run  procedures  outlined 
previously.  No  thickening  of  fin  surfaces  during  the  run  was 
observed.  The  color  of  the  R-113,  at  the  completion  of  the 
run,  was  amber.  The  bundle  data  (Figure  5.9)  does  demonstrate 
that  the  results  of  CNFT03  were  repeatable  at  least  within  the 
calculated  uncertainty. 

Single  tube  runs  were  made  in  accordance  with  the 
procedures  discussed  previously.  Nine  data  sets  were  taken  on 
each  tube,  at  the  same  conditions  specified  in  the  bundle  run. 
No  thickening  of  fin  surfaces  was  observed  during  the  runs, 
and  the  R-113  color  remained  amber.  There  is  no  major 
disparity  between  the  single  tube  results  (Figure  5.10),  and 
the  single  tube  results  for  CNFT03,  given  the  calculated 
uncertainty  band.  The  single  tube  results  were  significantly 
higher  than  bundle  performance,  and  the  build  up  of  non¬ 
condensable  gases  during  bundle  operation  and  evacuation  of 
the  apparatus  prior  to  commencing  single  tube  runs  is 
conjectured  as  a  possible  explanation.  At  no  time  during 
bundle  or  single  tube  operation  were  indications  of  non¬ 
condensable  gases  detected. 
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VEL.  (M/S) 

Figure  5.9 — CNFT04 — Tube  Bundle  Performance 
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VEL.  (M/S) 

Figure  5.10 — CNFT04 — Individual  Tube  Performance 


6 .  Smooth  Tubes  with  R-113  and  High  Coolant  Temperature 
SMT03  was  a  data  run  made  on  the  same  aforementioned 
copper  tubes,  but  with  R-113.  Vapor  saturation  temperature 
was  56.7  C.  Coolant  inlet  temperature  ranged  between  -1  C  and 
+6  C.  The  experiment  was  conducted  following  a  test  run  with 
R-113  and  finned  tubes  and  subsequent  replacement  of  the 
finned  tubes,  and  evacuation  to  11  in  Hg.  While  the  tube 
change  was  being  made,  the  R-113  remained  in  the  boiling 
chamber  exposed  to  the  atmosphere  for  a  two  hour  period.  The 
R-113  appeared  light  amber  in  color  at  the  commencement  of  the 
bundle  run.  The  desired  gage  pressure  was  maintained  by 
controlling  the  boiling  chamber  heaters  through  their 
corresponding  variac  controllers.  The  auxiliary  condenser  was 
not  utilized  because  the  additional  coolant  flow  provided  by 
the  auxiliary  condenser  made  it  difficult  to  maintain  the 
desired  system  pressure.  During  the  data  run,  the  top  two 
instrumented  tubes  assumed  an  orange  peel  textured  or  dimpled 
appearance.  The  R-113  had  a  darker  amber  appearance  at  the 
completion  of  the  bundle  run.  Again,  data  were  taken  at  the 
same  nine  coolant  flow  settings  following  the  procedure 
outlined  previously.  Data  taken  during  bundle  operation 
(Figure  5.11),  while  falling  within  the  acceptable  uncertainty 
bands,  clearly  demonstrates  the  presence  of  non-condensable 
gases  as  well  the  presence  of  the  contamination  acting  to 
inhibit  the  overall  heat-transfer  performance.  The  top  tube's 
very  poor  performance  can  be  attributable  to  its  location  in 
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Fiaure  5.11  —  BUT 03 — Tube  Bundle  Performance 


a  pocket  of  non-condensable  gases.  And  the  bundle  performance 
as  a  whole,  when  taken  in  comparison  with  the  performance 
demonstrated  in  the  bundle  run  of  SMT02 ,  indicates  an 
increased  concentration  of  the  contamination  and  the  presence 
of  non-condensable  gases. 

Upon  completion  of  the  bundle  run,  coolant  flow 
through  the  instrumented  tubes  was  shut  off  and  gage  pressure 
maintained  by  controlling  flow  through  the  auxiliary 
condenser.  While  no  coolant  flowed  through  the  instrumented 
tubes,  the  textured  appearance  noted  during  the  bundle 
operation  gradually  transformed  to  small  droplets  on  the  top 
two  tubes  which  then  coalesced  into  larger  ellipsoid  shaped 
droplets  that  were  held  on  the  tube  by  surface  tension  forces. 
Data  was  taken  on  the  first  two  tubes  in  the  bundle  following 
the  procedure  outlined  previously  for  single  tube  runs.  At 
the  completion  of  the  second  run,  it  was  noticed  that  the  top 
view  port  flanges  were  cold  to  the  touch  while  the  bottom  view 
port  flanges  were  warm  to  the  touch.  The  assumption  was 
therefore  made  that  non-condensable  gases  had  collected  in  the 
top  of  the  condensing  chamber.  Before  evacuation  of  the 
apparatus  could  be  effected,  the  pyrex  glass  in  the  top  center 
view  port  cracked.  As  a  result,  the  apparatus  was  secured  and 
no  further  data  taken.  It  should  be  noted  that,  prior  to  the 
view  port  breaking,  there  was  no  evidence  to  suggest  that  the 
apparatus  was  other  than  air  tight,  so  that  the  source  of  non¬ 
condensable  gases  was  most  probably  due  to  their  generation 
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within  the  boiling  chamber.  Data  taken  in  the  single  tube 
runs  (Figure  5.12)  reflects  not  only  the  effects  of  the 
contamination,  but  also  the  effects  of  the  presence  of  non¬ 
condensable  gases  as  well.  This  data  showed  the  greatest 
variance  between  calculated  and  observed  uncertainties  in  the 
overall  heat  transfer  coefficient.  This  data  is  limited  to 
the  performance  of  only  the  top  two  tubes,  due  to  the  afore¬ 
mentioned  breakage  of  the  view  port  and  the  presence  of  non¬ 
condensable  gases. 

C.  SUMMARY  OF  RESULTS 

The  comparison  between  finned  and  smooth  tubes  shows  an 
enhancement  ratio  in  the  overall  heat  transfer  coefficient  for 
the  finned  tubes  compared  to  the  smooth  tubes  of  approximately 
2.0.  While  this  enhancement  in  the  overall  heat-transfer 
coefficient  might  translate  to  a  vapor-side  heat-transfer 
coefficient  enhancement  ratio  of  approximately  4.0  to  4.5, 
based  upon  the  calculated  values  for  the  vapor-side  heat- 
transfer  coefficient  for  the  tested  smooth  tubes  when  compared 
against  the  Nusselt's  correlation  for  smooth  horizontal  tubes 
and  the  reported  enhancement  ratios  from  other  finned-tube 
investigations  [Refs.  5-13],  it  is  estimated  that  the  effect 
of  the  contamination  during  these  tests  has  been  a  degradation 
of  up  to  50%  in  the  vapor-side  heat-transfer  coefficient.  The 
contamination  has  resulted  in  the  promotion  of  another 
resistance  to  heat  transfer  that  has  no  predictable 
characteristics  and  can  not  be  quantified.  This  contamination 
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VEL.  (M/S) 

Fiaure  5.12 — SMT03 — Individual  Tube  Performance 


must  therefore  be  removed  before  any  quantifiable  heat- 
transfer  results  can  be  obtained. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Based  on  the  data  gathered  during  this  investigation  for 
condensation  of  refrigerants  R-114  and  R-113  in  the  multi-tube 
apparatus,  the  following  conclusions  are  reached: 

1.  The  construction  of  a  multi-tube  condensation  test 
apparatus  begun  by  Murphy  [Ref.  3]  and  Zebrowski  [Ref. 
4]  has  been  completed,  instrumented  and  integrated  with 
the  various  support  systems. 

2 .  Data  reduction  software  has  been  developed  that  produces 
high  confidence  results  from  experimental  readings. 

3.  Data  were  collected  from  condensation  experiments  with 
R-114  and  R-113  on  smooth  and  finned  tubes  demonstrating 
an  enhancement  ratio  of  approximately  2 . 0  in  the  overall 
heat  transfer  coefficient  for  the  finned  tubes  over  the 
smooth  tubes. 

4 .  Serious  apparatus  internal  contamination  has  degradated 
the  value  of  the  overall  heat-transfer  beyond  acceptable 
limits,  making  further  data  reduction  meaningless. 


B.  RECOMMENDATIONS 

Based  upon  the  data  gathered  during  the  operation  of  the 
multi-tube  condensation  test  apparatus,  and  the  experience 
gained  in  construction  of  this  apparatus,  the  following 
recommendations  are  made: 

1.  Samples  should  be  taken  from  the  R-114  and  R-113  used 
in  the  experimental  runs  and  a  chemical  analysis 
performed  on  these  samples  to  assist  in  determining  the 
source  of  the  contamination  in  the  system. 

2.  Modifications  to  the  apparatus,  as  dictated  by  the 
results  of  the  chemical  analysis  of  the  refrigerants, 
should  be  made. 
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3.  New  flow  meters  allowing  a  greater  range  in  coolant 
velocities  should  be  installed.  The  greater  range  in 
coolant  velocities  will  enhance  the  accuracy  of  the 
modified  Wilson  plot  method  in  calculating  the  outside 
heat-transfer  coefficient. 

4 .  The  temperature  control  system  for  the  eight  ton 
refrigeration  unit  should  be  re-evaluated  to  determine 
if  tighter  control  of  sump  temperature  can  be 
maintained. 

5.  After  appropriate  modifications  to  the  apparatus  have 
been  made,  the  smooth  tube  and  finned  tube  experiments 
should  be  repeated  and  contrasted  with  the  results 
presented  here,  in  order  that  baseline  tube  performance 
might  be  established. 

6.  During  operation  of  the  system,  flow  should  be  throttled 
to  the  coolant  supply  pumps  with  the  valves  on  the 
coolant  inlet  side.  This  will  increase  the  pump  life 
and  reduce  noise  in  the  laboratory  space. 
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APPENDIX  A 


UNCERTAINTY  ANALYSIS 


A  certain  amount  of  uncertainty  exists  in  any  engineering 
measurement.  These  uncertainties  arise  both  from  known 
sources,  such  as  calibration  and  measurement  errors  from 
sensing  devices,  and  also  unknown  sources,  such  as  operator 
experience  and  unexpected  experimental  reactions.  The 
procedure  used  in  this  uncertainty  analysis  is  based  on  the 
Kline-McClintock  [Ref-  21]  method.  This  method  assumes  that 
a  result  R  is  a  function  of  the  variables  that  contribute  to 
that  result.  Therefore  the  uncertainty  in  a  result  is  a 
function  of  the  uncertainties  in  each  of  the  variables.  This 
method  is  expressed  mathematically  by: 


3R  r,3R  .  3VV  .  /3R  .  3V2.2  i  1/2 
-R  "  ~] 1  +  (3vj  nr  + 


(A.l) 


The  measurement  of  the  overall  heat  transfer  coefficient 
(U0)  performed  in  this  thesis  is  given  by: 


TJ  =  _ — - 

o  A  -LMTD 
o 

where  the  heat  transferred  (Q)  is  given  by: 


(A. 2) 
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(A. 3) 
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and  the  Log  Mean  Temperature  Difference  (LMTD)  is  defined  as: 


IMTD  = 


T  -Tc . 

,  ,  sat  i, 

l0gtT  Tc  1 

sat  o 


(A. 4) 


Therefore  the  uncertainty  in  the  measurement  of  the  overall 
heat  transfer  coefficient  (U0)  is  given  by: 


=  !(«i)2  +  (^E)2  +  ,M)2  +  (!j2)2  +  (*^^,2,1/2 
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where — -  reduces  to — -  ,  and  the  uncertainty  in  the  LMTD  is 
A  D 

o  o 

given  by: 
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(A. 6) 


where  A  is  given  by  the  expression: 


AT  1  i 
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where  B  is  given  by  the  expression: 


(T  -Tc.)  T  , -Tc ■ 

sat  X  r  sat  ln 

log[T  -T~ 

sat  o 


where  C  is  given  by  the  expression: 
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It  should  be  noted  that  the  uncertainties  for  mass  flow 
rate  (r)  and  the  coolant  temperature  rise  across  the 
condensing  tubes  (At)  were  calculated  during  calibration  of 
the  flow  meters  and  thermopiles  respectfully.  Reference 
material  providing  data  on  thermophysical  properties, 
specifically  the  specific  heat  (Cp)  ,  listed  no  uncertainties 
associated  with  the  curves  for  aqueous-ethylene  glycol 
solutions.  The  composition  of  the  aqueous-ethylene  glycol 
solution  by  weight  percentages  was  calculated  by  measurement 
of  the  specific  gravity  and  the  change  in  this  composition 
over  the  experimental  period  was  determined  to  be  negligable 
and  therfore  this  uncertainty  has  been  ignored. 

The  uncertainty  in  the  overall  heat  transfer  coefficient 
represents  the  maximum  uncertainty  present  from  the 
calibration  of  the  four  flow  meters  and  the  maximum 
uncertainty  in  the  calibration  of  the  four  thermopiles, 
regardless  of  association,  and  therefore  represents  a 
conservative  estimate  of  the  uncertainty  present  in  each 
calculation  of  the  overall  heat  transfer  coefficient  (U0) . 
This  maximum  uncertainty  was  calculated  to  be  4.5  percent. 

The  disparity  between  the  calculated  uncertainties  and  the 
uncertainties  observed  by  the  spread  in  data  particularly  in 
the  single  tube  runs,  can  be  accounted  for  by  the  as-yet 
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undetermined  chemical  or  phase  reactions  generated  by  the 
contamination  in  the  system. 

The  values  of  the  uncertainties  in  the  measured  variables 
that  were  used  to  calculate  the  uncertainty  by  the  Kline- 
McClintock  method  are  listed  in  Table  A.l. 


VARIABLE 


TABLE  A.l 

UNCERTAINTY  VARIABLES 

VALUE  REMARKS 


6T 

r 


0.01  From  flow  meter 

calibration  data 


6  AT 

AT 


0.00  Not  available 


0.01  From  thermopile 

calibration  data 


oD 

o 


D 

o 


0.002  From  micrometer 

name  plate  data 


c  LMTD 
LMTD 

SMT02 

CNFT01 

CNFT02 

CNFT03 

CNFT04 

SMT03 


Listed  below  by 
run  title 

0.005 

0.002 

0.001 

0.003 

0.002 

0.017 
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APPENDIX  B 


LIGHTOFF  AND  SECURING  PROCEDURES 

A.  SYSTEM  LIGHTOFF 

1.  Push  the  starter  button  in  the  control  box  for  the  re¬ 
circulation  pump.  This  control  box  is  located  on  the  bulkhead 
above  the  re-circulation  pump  in  the  outside  area  adjacent  to 
the  refrigeration  unit. 

2.  Turn  the  switch  on  the  refrigeration  unit  control 
panel,  located  in  front  of  the  refrigeration  unit  to  the 
"auto"  position  after  passing  through  "on"  position. 

3.  Set  the  desired  temperature  on  the  roughly  graduated 
Fahrenheit  scale  on  the  control  panel  thermostat.  It  requires 
approximately  four  hours  to  chill  the  sump  40  degrees  C.  The 
thermometer  located  on  the  side  of  the  ethylene  glycol  sump 
must  be  monitored  to  ensure  the  desired  sump  temperature  is 
attained  and  maintained.  Slight  adjustments  in  the 
refrigeration  unit  thermostat  can  be  expected  due  to  the 
coarseness  of  its  scale. 

4.  Energize  the  heater  variacs  by  switching  on  the 
breakers  in  the  breaker  panel  located  in  the  laboratory  space 
on  the  bulkhead  next  to  the  counter. 

5.  Set  the  heater  variacs  to  the  desired  position,  after 
ensuring  that  the  switch  panel  for  the  heater  tubes  located  on 
the  apparatus  has  all  switches  in  the  "on"  position.  Monitor 
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apparatus  pressure  through  the  pressure  gage  at  the  top  of  the 
apparatus,  ensuring  system  pressure  does  not  exceed  30  psig. 

6.  Turn  on  the  pump  motors  by  pushing  down  on  the  arm  of 
the  appropriate  breaker  box  for  the  pumps  located  on  the 
bulkhead  next  to  the  ethylene  glycol  sump.  The  pumps  are 
marked  "auxiliary  condenser"  and  "instrumented  tube 
condenser,"  respectively.  Flow  in  the  auxiliary  condensate 
system  can  be  controlled  with  the  individual  gate  valves 
located  at  the  coil  penetrations  on  the  apparatus.  The 
auxiliary  condenser  will  produce  the  fastest  adjustments  to 
system  pressure  if  pressure  is  rising  too  quickly.  Flow 
through  the  instrumented  tubes  can  be  controlled  by  the  ball 
valves  located  at  the  bottom  of  the  respective  flow  meter. 

7.  Throttle  down  coolant  flow  to  the  supply  pumps  with 
the  valve  located  before  the  pump  suction.  This  will  increase 
the  pump  life  and  reduce  noise  in  the  laboratory  space. 

B.  SECURING  PROCEDURES 

1.  Turn  all  variacs  to  the  zero  position  and  switch  off 
all  breakers  in  the  power  panel  on  the  bulkhead. 

2.  Turn  the  breakers  for  the  pumps  to  the  off  position  at 
the  switch  boxes  near  the  ethylene  glycol  sump. 

3.  If  apparatus  will  not  be  operated  for  an  extended 
period,  turn  the  switch  on  the  refrigeration  control  panel  to 
the  "off"  position  after  passing  through  "on." 

4.  Allow  the  re-circulation  pump  to  operate  for  at  least 
ten  minutes  after  switching  off  the  refrigeration  unit  to 


73 


APPENDIX  C 


1000  ' 
1005  ' 
1010 
1015  1 
10Z0 
1025  1 
1030 
1035 
1040 
1045 
i050 
1055 
1050 
1065 

1070 

1075 

1090 

1065 

1030 

1095 

1100 

lies 

1110 

1115 

1120 

1125 

1130 

1135 

1140 
1145 
1150 
1 155 
1160 
1165 
1170 
1 175 
1180 
1185 
1190 
1195 
1200 
1205 
1210 
1215 
1220 
1225 
1230 
1235 
1240 


DATA  REDUCTION  PROGRAM 


FILE  DRP1F 

cURP0SE  This  program  collects  and  processes  condensat lcn  data 


the  R-114  tuOe-pundle  apparatus. 
CREATED  N0UEM8ER  2.  1988 
UPDATED  Novemper  6,  1986 
CHANGE  H00=3  ANO  USING  CALCULATED  Cl  FROM  WILSON 
BEEF 

PRINTER  IS  1 


PRINT  USING  "4X . " 'SELECT  OPTION 


PRINT  USING  ''EX,'  “ 0 
RRINT  USING  "  6  '  ,  '  1 
FPIN7  USING  '  6 X  .  2 

PRINT  USIN6  "EX, -''3 
PRIN*  USING  'GX,— 4 
PRINT  USING  " 6 X  ,  " " 5 
PRINT  USING  “EX,"”S 
PRINTER  IS  701 


TAKING  DATA  OR  REPROCESING  PREVIOUS 
PLOTTING  h  US  DELTA-T' ' ' 

PLOTTING  HRAT  US  N""“ 

PLOTTING  WILSON"'" 

PURGE  RILES"'"' 

XVREAD . 

NuSSE_"  ESTIMATE . 


INPU*  Icall 


IF  Icall-0  THEN  2Aul  Main 

IF  leal  1  =  1  THEN  CALL  P i 0 1 2 

IF  Icall “2  THEN  CA_L  Pletl 

IF  leal i*3  THEN  CAuL  P i c 1 3 

IF  Icall“4  THEN  CAuL  Purge 

IF  Icall=5  THEN  CALL  Xyreac 

IF  Ica!l*6  THEN  CA^L  Nusseit 


CAT  f. 


END 


SUB  Main 
COM  /Cc/  C( 7  > 

COM  /Fid/  I  ft 

COM  /Nus/  lin.Tsat .Odd ,Hnus 

COM  /Will/  Doa' 4  ) , Di a f 4  / , Kma ■ 4  ) , lac t 

COM  /Wi 1 2/  Delta, Isat.Nsets, Hod, Cia(3),Aipaai3> 

DIM  Fma(8,4l,Fm(4),Emf(6,Tpi3i,T(8),Ho(3),0dp( 3),Uo<3) 

OATA  33.4.20,34.4,34.5,23.5 

DATA  41.3,25,41.3,41.6,28.8 

DATA  49.1,30,48.4,48.6,34.1 

DATA  57.0,35,55.3,55.6,38.4 

DATA  64.8,40,62.3,62.7,44.6 

DATA  72.7,45,69.2,69.7,48.8 

OATA  80.6,50,76.3,76.7,55.2 

DATA  88.4,55,83.3,83.8,60.5 

DATA  96.3,60,90.2,90.8,65.6 

READ  Fma  < •  ) 

OATA  5.172,5.172,5.172,5.172 
READ  C 1 a< •  ) 

DATA  0.10086091,25727.94369,-767345.8295,78025595.81 
DATA  -9247486589,6. 976 8 8E 1 1 , -2 . 66 1 92E13 , 3 . 94078E1 4 
READ  C< •  ) 

DATA  0.015875,0.0 14000 ,0.0, 0.0, 0.0 
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;  _4;  un  i  H  0  .  01  0LS3  ,  0  .  U:  0ib0  ,  0  .  0  .  0  .  0  .  0 .  t“ 

1250  DATA  336. 0, 42. S^S, 0.0, 0.0, 0.0 

1255  DATA  B.00055S9.3 

i  DEC  READ  Doa 1  •  > . Dial  ‘  i  ,Kna>  »  '.Delta, hoc  1  -c  d=n  •'  C : 

12E5  L-1.21S2  1  Condensing  length 

1270  Jset-O 

1275  BEER 

1280  INPUT  “ENTER  MONTH,  CATE  AND  TIME  (MM  CD  Hm  •  MM  S3  ' “ , C  t  gl 
1255  OUTPUT  709. "7C" ;DtgS 
1290  BEEF 

1295  INPUT  “SELECT  OPTION  (0-OAQ,  1 “FILE  )  “  , In 

1300  Iha^c*! 


1305 

1310 

1315 


BEEP 
INPUT 
BEEF 
INPUT 
I  ip«  1 


"WANT  A  HARDCDFV  PRINTOUT  1  =DEF-'r  ES  ,  0  =  N  >  *  ,  I  hard 
“SElECT  (0-R-H4. 1=STEAM,2=p-1I3,3*EG  )"  .  I  ft 


1330 


I sat  =  2 


1340 

1345 

1350 

1 3S0 
1 3E5 
1 370 

IC^E 


BEEC 

INPUT  “SELECT  SAT  TEMP  MODE  < 0=L ID . I =VAp . 2  =  ' L 1 0--VAF l.'D-DEF  i " , I  sat 
:c  Ina-3-i  Then  PRINTER  IS  701 
IF  >  =  0  THEN 
BEEP 

INPUT  “GIVE  A  NAME  PGR  THE  NEW  DATA  FIuE", PiieS 
CREATE  BOAT  File*, 20 
ASSIGN  @F i le  TO  File* 

BEEC 


INPUT  “ENTER  TUBE  CODE'.Itube 
BEEP 

INPUT  "ENTER  EG  CONCENTRATION  •.  WT  PEPCENT  ) “ , Egrat 
ENTER  709 , D t gS 
OUTPUT  @F i le . Otg$ 

OUTPUT  @p lie-  I  tube  ,  Egret .  Ddi . DdC . Dd3 . Dd4  .  Dd5 
Iact-10 
BEEF 

INPUT  '  SEl.EC 7  ( 0  =  TCF , 1 -SECOND . 10=5UNDlE=DEFAuLT : " ,  last 

PRINT 

PCINT  USING  “ 1 0 X , “ " c ! L E  NAME  '  .ICA'.File* 

P  R I N  T 
1440  ElSE 
1445  BEEF 

1450  INPUT  “ENTER  NAME  OF  EXISTING  FILE', File* 

1455  BEEP 

14E0  INPUT  “ENTER  NUMBER  OP  DATA  SETS  STORED “, Nset s 

14E5  I w i 1=1 

1470  BEEP 

1475  INPUT  “WANT  TO  CALL  WILSON7  < 1 -DEFAULT- YES . 0=NO > " , I u l 1 

1480  IF  Iwll-1  THEN 

1485  8EEP 

1490  INPUT  "WHICH  TUBE  ( 0=TOP , 1 -SECOND ....  10-BUNDLE  ) "  ,  lac t 

1495  CALL  Wilson 

1500  IF  I  hard- 1  THEN  PRINTER  IS  701 

1505  PRINT 

1510  PRINT  USING  "10X,""FILE  NAME  " " , 12A" . Pi leS 

1515  PRINT 

1520  IF  lact  >9  THEN 

1525  PRINT  USING  “ 10X , " " INSIDE  COEFFICIENTS-  " “ , 4 ( DD . 3D , 2X > “ . C i a > • > 

1530  PRINT  USING  ' 1 0X ALPHAS  " ” , 4(1 X , M2 . 3DE , 2 x  '.A’D 


1 3S5 
1390 
1395 

1400 

1405 

1410 

1415 

»  *  -  n> 

:  *zs 

1430 

1435 
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iszs  else 

1540  PRINT  USING  " 10X , ' ‘ INSIDE  COEFFICIENT  FOR  TUBE  "‘.D, 

D '  . lac t , C; a< lact  ) 

1545  PRINT  USING  ” 10X , " "ALPHA  FOR  TUBE  " "  ,  D , " "  ; 

C" . lact , Alpaai lact  ) 

1550  END  IF 

1555  ELSE 

1565  INPUT  "ENTER  Cl  VALUES  ( DEFAULT =5 . 172  > " , Dial ♦  ) 

1570  END  IF 

1575  ASSIGN  9File  TO  File* 

1580  ENTER  @FiIe,Dtg$ 

1555  ENTER  @F i le . I t ube , EQra t , Ddl ,Dd2 , Dd3 . Dd4 , Dd5 

1586  INPUT  "ENTER  TUBE  CODE",Itube 
1 5S0  END  IF 

1595  lout-1 

1600  BEEP 

1605  INPUT  " WANT  TO  CREATE  AN  OUTPUT  FILE7  < 1«DEF=YES , 0-N0 >"  lout 
1610  IF  I ou  t  *  1  THEN 

1615  BEEP 

1620  INPUT  "ENTER  A  NAME  FOR  OUTPUT  FILE'.FoutS 

1 E 0 5  CREATE  BOAT  Pouts, 5 

16.10  ASSIGN  @Fout  TO  FoutS 

1 E 35  END  IF 

164  0  i 


1546 

Do-Deat I  tube  > 

1  ESC 

Di"Diaf I tube  > 

1555 

Km*kna< I  tube  ) 

1660 

A*»PI*Di "2/4 

1  Cross-sectional  area 

1665 

Ao”PI«Do*L 

1670 

Rr»«Do»LOG(Do/Di  )  /  (  2  m  ) 

1E79  i 

1660 

IF  Im=0  THEN 

1555 

PRINTER  IS 

1 

1690 

BEEP 

1655 

PRINT  "SET 

FlOUME'E=  READINGS  CORRESPONDING  TO" 

1~00 

PRINT  " 

" .Fna: Jset , 1); "%  OF  METER  2  AND  hIT  CONTINUE 

1705 

PAUSE 

:-i0 

PRINTER  IS 

701 

1 "  1 5 

OUTPUT  709. 

"AR  AP0  AlB  UR5 ' 

1-20 

Nerd=3  i 

INCREASE  TO  9  IF  FIVE  TUBES  IN  BUNDLE 

2  -  -  c 

FOR  1=0  TO 

Nend 

OUTPUT 

709, ’AS  SA¬ 

1725 

Vsun=0 

1740 

FOR  J*1 

TO  5 

1745 

ENTER  709. E 

1750 

Vsur'*V3ufy'+E 

1 755 

NEXT  J 

1760 

Enf  (  I  >» 

‘Vsum/B 

1765 

NEXT  I 

1  770 

OUTPUT  709, 

“AR  AF20  AL23  UR5 " 

1775 

FOR  1-0  TO 

3 

1780 

OUTPUT 

709; "AS  SA¬ 

1785 

VsuM-0 

1790 

FOR  J-l 

TO  50 

1795 

ENTER  709, E 

1800 

Usun-Usun+E 

1805  ' 

WAIT  .25 

1610 

NEXT  J 

16:5 

Tp< I  )=Vsum/ 50 

«  r*  *h  -> 
iObi1 

NEXT  I 

1825 

ELSE 

1830 

ENTER  @Fi 2  e . Fm( »  )  ,Enf ( »  > ,  Tp  t  »  > 

1835 

END  IF 

1840  ' 

1845  ' 

DATA  ANALYSIS 

1850  ‘ 

1 8S5 

Nend=S 

1860 

FOR  1*0  TO  Nend 

1865 

T( I )“FNTvsv(Emf(  I  ) ) 

1870 

NEXT  I 

1875 

Tvap*<  T (  0  )  +  T (  1  )+T(  2  .}  1/3 

1880 

.  T1  iq*( T< 3 >  +  T (4)1/2 

1890 

Tsat »T1  iq 

1895 

END  IF 

1896 

IF  Isat-1  THEN 

1900 

T  sa t  =  Tvap 

1905 

END  IF 

1905 

IF  Isat-2  THEN 

1907 

Tsa t «( Tvac*Tl ; -  i/2 

1906 

END  IF 

1910 

FOR  1=0  TO  4 

1915 

Fwt  I  i=Fmaf  Jset  ,  I  • 

1920 

NEXT  I 

1925 

Jset-Jset+1 

1920 

PRINT- 

1925 

PRINT  USING  "10X,""Oata  set  njmper  -  ”  “  ,  DD "  . J se t 

1940 

PRINT 

1945 

lbeg=0 

1950 

lenci  =  3 

1955 

IF  Iact  -10  THEN 

I960 

Ibeg= Iact 

1955 

lend* I ac  t 

15"0 

END  IF 

1  9"5 

FOR  I*Ibeg  TO  Ieno 

1990 

Gnad  =  FNG'_ad(  Emf  (  I  =5  ;  1 

1  965 

De  1 1  *A8Si  Tp  (  I  '  Grad*  10.  ) 

1990 

Tavg=T< 1+5  ) - D e 1 1 ♦ .5 

1995 

Rhoeg*FNRhoeg' Tavg,Egret ) 

2000 

Nueg=FNNueg) Tavg,Egrat  > 

2005 

Mueg*Nueg»Phoeg 

2010 

Cpeg*FNCpeg(Tavg,Egnat  :■ 

2015 

Keg*FNKeg( Tavg , Egrat  ) 

2020 

Pneg=Cpeg*Hueg < Keg 

2025 

Hdot  *FNFmca j  < I , T( 1*5 1 ,  rmi I  ;  ) 

2030 

Ueg*Mdot /( Rhoeg*A.* > 

2035 

Reeg=Ueg*Oi /Nueg 

2040 

Res*4«Mdot/'PI«Mueg»' 0i-4  • Delta  )  ) 

2045 

Qdot=Mdot»Cpeg»Delt 

205C 

Qdp( I  )=Qdot / Ao 

2055 

IF  1*0  OR  I=Iact  THEN 

2060 

Qdp 1 =0dp ( I  ) 

2065 

CALL  Nusselt 

2070 

END  IF 

2075 

Lmt d*Oe 1 1 /LOG< ( Tsat -T ( I +5  )  )/ ( Tsat-T( 1+5  >-De 1 1  )  ) 

2080 

Uo<  I  )*Qdp  (  I  )  /Lnt  d 

2065 

IF  Reeg<  4000  THEN 

78 


.Cat 

2095 

2100 

2:0s 

:::c 

Z1 15 


Nueg  =  C  is'  I >*t 1*5. 434E-3+preg'  .?•  <  Res/  Hod  )'  1 .25  >‘  .5 

ElSE 

BEEP 

PRINT  USING  “10X, " -INCORRECT  TURBULENT  CORRELATION" 
Nueg= . 027+Reeg’  .8+Preg'  .3333*Cfeg 
END  IF 


2120 

Hi 

«Nueg*Keg  xQ : 

2125 

Hot  I  >  =  1  /  t  1  /UCH  I 

>-DcX 

Di»H:  >-Rm  ) 

2130 

IF 

1=0  OR 

I » I act  THEN 

2135 

PRINT 

USING 

“10X, 

"Mass  flow  rate  * 

“ " .MZ.3DE" 

,  Mdot 

2136 

PRINT 

USING 

"10X, 

'"Inside  Tube  Dia.  (m.)  ” 

"“.MZ.3DE" 

,Dia( Itube 

; 

PRINT 

USING 

“10X, 

'"180  DEG  OUER  Dia.  ( HOD  )  - 

" " ,MZ .  3DE" 

;  Hod 

2140 

PRINT 

USING 

"10X, 

'"'Inlet  temperature 

"" .MZ.3DE" 

;T< 1  +  5  ) 

2145 

PRINT 

USING 

"10X, 

"'Saturation  temp  (Deg  C>  • 

■ " , MZ . 3DE " 

.  Tsat 

2  1  46 

PRINT 

USING 

■■  10X  , 

""DELTA  Tape  Thickness  * 

" " ,MZ .3DE" 

, Delta 

2147 

PRINT 

USING 

"10X. 

"  "DELT  ten p  Di  f .  “ 

"\MZ.3DE" 

;  Del  t 

2148 

PRINT 

USING 

n10X, 

““Log.  Mean  tenp  Dif. 

” " ,MZ .3DE' 

.  Lmt  d 

2150 

PRINT 

USING 

“10X, 

"“Heat  flux  “ 

"" .MZ.3DE' 

. Qdp  t  I  ) 

2151 

PRINT 

USING 

"10X  , 

""Conductivity  E.G.  " 

""  .MZ.3DE" 

.  Keg 

-  ■  C£ 

PRINT 

USING 

-  10X  . 

""Conductivity  Tube  Metal  ■ 

“ " , MZ . 3DE ' 

;  K  m 

2155 

PRINT 

USING 

"10X, 

" "Frandt 1  number  = 

' “ , MZ . 32E ' 

,  Preg 

2 1  £0 

PRINT 

USING 

"10X, 

""Reynolds  nunper  = 

'  MZ.3DE' 

.  Reeg 

2151 

PRINT 

USING 

"10X. 

""Reynolds  number  H&B  S  * 

" " , MZ . 3DE " 

,  Res 

2165 

PRINT 

US  I NG 

"10X. 

""Inside  h.t.c.  “ 

" “ , MZ . 2DE ' 

■  Hi 

2166 

PRINT 

USING 

"  10X  , 

""Inside  NUSSULT  NO.  » 

" " , MZ . 3DE ' 

.  Nueg 

2157 

PRINT 

USING 

••  10x  , 

""OUERALL  H.t.c.  <  Uo  > 

" “ ,MZ .30E' 

.  Up  ■  I 

2170 

PRINT 

2  175 

PRINT 

USING 

"  10X  , 

""TUBE  FM  UEG 

DELT  Uo 

hO 

HNUS 

..  « 

2180 

PRINT 

USING 

"10x , 

""  *  ( X  x  tn/si 

(  K  ) 

(  Ui  /  cm  ’  2  . 

K  ; . 

2185 

PRINT 

USING 

"  10X  , 

30 , 4X , 30 . DO . 3X . Z . DO , 4X , DO . 30 , 2X , Ml . 3DE . 3  X , m . 32  E 

.3* ,m: 
Z1SS 

7  1  Q- 


,11 

:cs 


30E  ,  3  * " .  I  *  1 ,  Ft*'  I  '• ,  Ueg  .  De  1 : .  Uc '  I  ) ,  He '  I  > ,  Hnu  s 

else 

PRINT  USING  "  10>  ,  3D  ,  4x  ,  32 . 02 . 3>  ,  Z  .  DD  .  4X  ,  DD  .  3D ,  2>. ,  MZ  .  325 , 3>  ,MZ  .305 
►  1 , Pm (  I  ) ,  Deg ,  De  1 1 ,  Jo '  I  ; .  He  •  I  • 

ENO  IF 
NEXT  I 

IF  >«e  Then 


ZZS5 

2290 

2295 


Ox  accpt=l 
BEE= 

INpUT  "OX’  TO  ACCECT  TUI3  cg-r  (  i -DEFAUlT-TES  ,  0  =  NO  P "  ,  Ox  a :  cc  t 
IF  Qxaccp".*!  THEN  OUTPUT  @F  i  le  .  Fmi  •  ) ,  Emf  '■  *  ) ,  Tpt  •  ) 

END  IF 

IF  (Okaccpt=l  OR  Im=i)  AND  Iout«l  THEN 
FOF  1=0  TO  3 

OUTPUT  @Fput.HO‘I  ,0dpi  I  > 

NEXT  I 
END  IF 

IF  Ih»0  THEN 
OX  rp  t = 1 
BEEF 

INPUT  "DILL  THERE  BE  ANOTHER  RUN  <  1  =  VES*=0EF AULT  ,  0-NO  > ’  ,  Ox  rp \ 
IF  OX  rp t  =  1  THEN  1660 

Else 

IF  Jset-'Nsets  THEN  1680 
END  IF 


2305  ASSIGN  @F i I e  T0  • 
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2710  IF  Iout=l  THEN  ASSIGN  @Fiie  TO  * 

2215  SuBENO 

2320  DEF  FNG-ad( T  ) 

2325  Graa  =  -3.877857E-5-2*4.7l4285':E-S»T 

2330  RETURN  G^ad 

2235  FNEND 

2340  DEF  FNKcu< T  ) 

2345  •  OFHC  COPPER  250  TO  300  K 
2350  Tk  -Tf-273 . 15 

2355  K  =  434- . 1 1 2  *  T  k 

2360  RETURN  K 

2365  FNEND 

2370  DEF  FNNueg ( Tc , Eg^ > 

2375  i  RANGE  OF  UALIDITY  =  -20  TO  20  DEG  C 

2380  Tk-Tc+273.15 

2385  Nul*7.1196507E-3-Tk»(7.4863347E-5-Tk»<2.6294943E-7-Tk *3 . 0833329E- 1 0  ' 

2390  Nu2*4.9237E38E-3-TM<4.9213912E-5-Tk«< 1 . 6437S34E-7-Tk »1 . 8333331 E- 10 ) 

2395  Nu3-8.E58E293E-3-Tk*(8.8837902E-5-Tk*<3.0495032E-7-Tk»3.4999996E-10) 

2400  A2-!Nu3-2«Nu2+Nul >/200 

2405  Al-< Nu2-Nul-940»A2  )/10 

2410  A0«Nul-42*Al-17E4*A2 

2415  Nu-A0  +  Egr» < A 1 +Egr« A2  ) 

2420  RETURN  Nu 

2425  FNEND 

2430  DEF  FNCpeg< Tc .Egr  ) 

2435  i  RANGE  OF  UALIDITY ■  0  TO  20  DEG  C 

2440  Tk-Tc+273.15 

up  i*i  .  «a /w isb(flt  +  i+ i  k»e  .  j 
2450  Cp2*l . 47461 25E*3+Tk  «G . 25 

2455  Cp3*9.S800500E+2+Tk*7.3 

2460  A2-( Cp3-2«Cp2+Cpl  >/200 

2465  Al«(Cp2-Cp 1 -900+A2  )/ 10 

2470  A0-Cpl-40.A1-1G00«A2 

2475  Cp-A0+Egr» i A 1 +  Egr» A2  ) 

2460  RETURN  Cp 

2485  FNEND 

2490  DEF  FNRhoeg !  *  .  Eg1"  ) 

2455  Pol-1 .0607093E+3-T«(3.70212E3E-i+T«4.0S3-'182E-2 

2500  Rc2=i .0748272Ei3-T*( 4 , 4266 1 95E-1 »4 . 0S3S705E-2 

2505  Po3=l .0885934E+3-T«i5.73S5653E-i*T«6.12£1405E-3 

2510  A2*(Ro3-2*Ro2+Rol)-'20O 

2515  A1 = ! RoZ-Rol -900+A2 > ■  1 0 

2520  A0-Rcl-40«A1-1600*A2 

2525  Ro«A0»Egr» ! Al+Egr»A2  ) 

2530  RETURN  Ro 

2535  FNEND 

2540  OEF  FNPreg'  T  ,  Egr  s 

2545  Pr-FNCpeg‘  T.Egr  )+FNNueg!T,Egr  >«FNRceg  T  .  Eg1-  1  /  FNk'eg  ■  T  ,  Eg1-  ■ 

2550  RETURN  Pr 

2555  FNEND 

2550  DEF  FNKeg' Tc , Egr  ) 

2565  i  RANGE  OF  UALIDITY  -20  TO  20  DEG  C 

2570  Tk-Tc+273.15 

25"5  Kl-2 .2824709E-1  +  T1 •( 5.5989285E-4  +  Tk  *3 . 57 1 428EE-7  > 

2580  K2»2.5846616E-1+Tk*(2 . 3978571E-4+T1 +7. 1428571E-? ; 

2585  K3-3 . 2 138932E-1 -Tk • i 3 . 0042857E-4-Tk • 1 . 428571 4E-5  i 

2590  A2-!  K3-2  +  k2+k:  1  5/200 

2595  Al-( K2-K1-900.A2  )/10 

2500  A0-K1-40.A1-1S00+A2 

2505  K • A0  +  5gr » (  a  1 *Egr • A2  : 
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:b:c 

6  E  -  L1 6  \  ► 

ISIS 

FNEND 

3B30 

DEr  ■  ' 

3  B  Z  E 

=  =£>  P ■  » 

ZB  30 

0-1  -P 

Z63E 

T  a  n  h  -  ■  p  -  Q  ■■  ■  t  P  -1  G 

2640 

RETURN  Tanr 

2645 

FNEND 

2650 

DEF  i 

2B55 

COM  /Cc  •'  C  '  "  ) 

26S0 

T  =  C ( 0  > 

2665 

FOR  1=1  TO  - 

2670 

t-t*c:  I  •*2"  I 

2675 

NEXT  I 

2680 

RETURN  T 

2685 

FNEND 

2890 

DEF  FNBe t  a ( T  ) 

2695 

Rop=FNRhc<  T  +  .  1 

2700 

Rom=FNRho ( T-  .  1  > 

2705 

Bet  a  =  -2  7 '  Pop  •‘■Rom  ••'Rop-Rom  it 

'  ? 

Z"7 1 0 

RETURN  Bets 

2715 

FNEND 

2720 

OEF  FNPsa  t i T c  ) 

2725  ' 

0  TO  80  deg  F  CURUE  FIT  OF  Psst 

2730 

Tf**l  . 8  *  T  c  m  32 

2735 

Pa=5. 945525-Tf ♦( . 153S2082+T f • < 1 .4840963E- 

2740 

Pg-Pa-1 4 . 7 

2745 

IF  P  g  >  0  THEN  i  +-PSIG,- 

-•in  Hg 

2750 

Psat  =Pg 

2755 

EL5E 

2760 

Psa t =pg»29 . 92/14.7 

27ES 

END  IF 

2770 

RETURN  F  s a t 

^  1  I'S 

r  lid'll.' 

2790 

OEF  FNFmcal'I.T.Fm) 

2785 

DIM  B 1 <  4  , E2\ 4  ;,B3i4  ),M1(4) 

,M2  ' 4  >,M3( 4  ; 

2-90 

DATA  2.48835E-2.6.71749E-3, 

-5.6103E-3.-5 

2795 

DAT 4  -i  .4'7621E-2,2.5320_’E-3 

, -2 . 75396E-3 ,  • 

2796 

DATA  1.048E-3.1.012S2E-2.1. 

1 14E-2 . -3 . 867 

2800 

DATA  1 .71525E-3.2.70428E-3. 

1.937496-3,1, 

2805 

DATA  2 .034E2E-3.2.86716E-3. 

2.09S17E-Z.Z. 

2806 

DATA  1 .^SeiE-S.Z.TeEZZE-Z.Z 

.  0012  IE-3 .2 . 2i 

2610 

READ  Bl'  •  '  ,B2<  •  >,B3'  •  >,M1 !  • 

) ,  M2  *  ) , M3 ' •  ) 

2811 

IF  T •:  -7  THEN 

2615 

Zl-61 '  I  )  +  Ml  (  I  )*F r, 

2520 

Z2  =  E2 <  I  >  +  M2(  I  >*Fm 

2621 

S*( Z2-Z1  )/ 13 

2822 

C  =  Z>8*20 

2824 

EL5E 

2825 

Zl-BZ'  I  .'-"M2!  I  )  »Fm 

2825 

Z2-B3'  I  )+M3(  I  /*Fi" 

262- 

S-< Z2-Z1  >/ 16 

2826 

C=Z if S*7 

2630 

END  IF 

2835 

Mdot -C  +  9*  T 

2840 

RETURN  Mdot 

2845 

FNEND 

2850 

SUB  Xy-ead 

2855 

BEEF 

Z-T f  •  9 . E 150E71E-6  )  ) 


.5079E-2.4.12SE-4 
-4 . 30913E-3 , -3. 10937 
17E-3.0.0 

92Z37E-3.Z .5S646E-3 
12433E-3.2.93127E-2 
240E-3.0.0 
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-if  '  f 

3075 
3030 
30S5 
3090 
3095 
3100 
3105 
3110 
3115 
3130 
3 1 35 
3130 
3135 
3140 
3145 
3150 
3155 


INP'JT  ’ENTER  FILE  NAME  "  ,  F  ;  I eS 
SEE- 

INPUT  "ENTER  NUMEEP  OF  X,Y  PAIRS", N 
ASSIGN  @c  1 1  e  tq  FiieS 
FOP  1  =  1  T 0  N 

ENTER  0F ; i e , x , V 
PRINT  X , V 

NE v  t  ; 

SUEEN2 
SUE  °-=ge 
BEEF 

I NPUT  "ENTER  FILE  NAME  TO  BE  DELETED ", F i 1 eS 

PURSE  FiieS 

GCTC  3910 

SUBEND 

SUB  Wilson 

COM  /Will/  Doe' 4  ) ,  D:av  4  > ,  Kr.a<  4  ) ,  lac  t 

COM  /Wil2/  De  1 1  a  ,  I  sa  t  ,  Nse t  s  , Hod,  C l  a i  3  ) ,  A]  paa i  3 

DIM  Fm 4 ) .Enf ( 5  i  ,Tp< 3  >,T< 8  >,Xa' 30  ) ,  Ya(20  ) 

BEEP 

INDUT  "PLEASE  RE-ENTER  NAME  CF  FILE", FiieS 
ASSIGN  @c l 1 e  TO  cile$ 

INFU”  "ENTER  TUBE  CODE", I  tube 
BEEP 

INPUt  "C-IUE  A  NAME  FOP  XY  FILE",Xy$ 

CREATE  BOAT  a  y  $  ,  5 
ASSIGN  @Xy  TO  X  y  $ 

L-i .3192 
Do«Doa' I  tube  > 

Di 'Dia ( I  tube ) 

Kn'Kna t 1 1 ube  > 

A»“PI*0i "3/4  1  Cross-sectional  area 

Ao«p I ♦Co*L 

Rr»*Do*LOG'Do/Di  )/<2«Km) 

Initial  values 
T  f  «  T  5  a  t 
Aiwa* . S55 
2  i  =5 . 172 
£'3.6! 
iceg-U 

I e n c  =  3  'CHANGE  TC  4.  IP  FIVE  TUBES  IN  BUNDLE 
IF  1  act  ,10  Tu£i\i 
I  beg* I ac t 
I ena* I ac  t 
END  IP 

FOP  I'lbrg  TO  lend 
S*  '0 
Sy-0 
S  *  s  »  0 
S*  y*0 
Jset'0 

ASSIGN  @p i 1 e  TO  FiieS 

ENTtF  @File.DtgS,Itube,Egrat,Ddl,Dd2,Dd3.0d4,Dd5 
ENTEC  @F 1 1 e ■ Fm •  ) , Enf 1 •  ) , Tp ( •  ) 

FDR  J=0  TO  6 

T ( J  )”rNTvsv( Enf  ! J  )  ) 

NE  X  T  J 

Tvap  =  i"’;0)  +  T'l  ;  +  T  (  2  !  5  ■'  3 
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3  1  9? 
3195 
3200 
3305 


3310  > 


334  0 
3345  ' 
3350 


3355 


3370 

CC^S 

3390 

3385 

3390 

3395  ■ 

3300  ' 

3305 

33  i  0 

3315 

3320 

3325  ' 

3330 


1  i  7^ 


3340 

3345 

3350 

3355 

3350 


3380 
33e5 
3395 
3400 
3405 
3410 
3415 
3430 
3435 
3430 
3435 
3440 
344  = 
3450 
3455 


’3;a-  T  3  '*7,  4  )  ;.’3 
IF  lsat*0  THEN 
Tsat-Tlig 
£■_  EE 

Tset  *Tvap 
END  IF 

6rad  =  FNG'-ao<  T (  1+5  )  ) 

Del  t  =  AES  ■:  Tp  ■  I  j/  (  Grao*  10  >  ) 

"avg*T i 1  +  5  1  +  De  i  t  *  .  5 

Ua  t  e+'/E  t  hy  i  ene  Glycol  Mixture  Properties 
Rhoeg-FNRhoeg1 Tavg.Egrat  ) 

Nueg*FNNueg>  Tavg.Egrat  , 

Mueg*Nueg*Rhoeg 
Cpeg=F N 3 peg< Tavg.Egrat  ) 
k eg=FNkeg ' Tavg.Egrat  ) 

Preg*Cpeg*Mueg/Keg 

Mdct*FNFncal(I,T;I+5,FniI  ;  ; 

0eg*Mdot / ( Rhoeg*A>  ; 

Reeg=Geg*0i .‘Nueg 

Pes“4»Mpot/( FI*Mueg*  Ei-4*Delta  .  ) 

Gdot*Mdot*Cpeg*0elt 

3dp*Qdot /mo 

Lntd-Delt/LOG; ( Tset -T , 1+5  >  )/ <  Tsa t -T < 1+5  )-De 1 
Uo*Qdc /Lnt  d 

Onega* (  l+S.484E-3*PregA  .7»\Res/Hoo  '  1 .25  .5 

R-114  Properties 
Hf g-FNHf q'.  T sat  > 

K  f "FNK  <  T  f > 

Rhof “FNRho1 T f  ) 

Mu  f *FNMu ; Tf  ) 

F-<Kf '3»Rhcf"2»G»Hfg/(Muf »Do  +  0dp  )  1*. 33333 

Hc*Alpa*F 

Two*T  sat-Odp/Ho 

Tf-Tsat/3+2*Tuo/3 

Y-<  1 /Uc-Rm  )*p 

x»0o*F'(Keg*0nega  > 

PRINT  *  OMEGA-"  .Onega  .  1  F-  "  F  ,  "x-  .  x  ;  "  y«"  :  Y 

Xa: Jset  >*X  '  INEFFICIENT  (MODIFY  LATER) 

Ya<  J set  )— Y 

Sx  *5*  +  X 

Sy-Sy+Y 

SxS=SxS+X«X 

Jset* j set +2 

IP  Jset'Nsets  THEN  3135 
ASSIGN  @p 1 1 e  TO  * 

S 1  ope* 1 Nse t s ’ S x y- Sx • Sy  )/ ( Nset s *Sx.s-Sx  '2  ’ 

I nt cp t - ( Sy-S 1  ope* Sx  )/Nsets 
Cic*l/Slope 
Alpac-l/Intcpt 
Cerr-ABSt ! C1-C1 c  ;/Cic  > 

Aerr-ABS; <Aipac-A]pa //Alpac ) 

IF  Cerr/,001  OR  Aerr.-.00l  THEN 
Aloa*<  Alpa  +  Alpac  > • .5 
C 1  *<  C 1 +C ic  )  *  .  5 

PRINT  ”  3 1 C- “ • C 1 c ; 1 ALPA* " . 1  pa 
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34EC 

2455 

34-0 

3i~c 

**  “  f 

3455 
3490 
34  95 
3500 
3505 
35:0 
35:5 
3523 

7C'  o 

3530 

3540 

3545 

2550 

3555 

3550 

35E5 


2595 

3590 

3595 

3800 

3605 

3810 

3815 

3620 

3625 

3630 

3535 

3640 

3645 


GC'2  2100 
EN2  :c 
EEEe 
EEEr 

FFIN'EF  15  1 

P=INT  :  C  1  c  ;  " ALFA* "  •  ft  1  pa 

2^a^  I  >=C  1  c 

Alpaai  I  '-Alpac 

PRINTER  IS  701 

FOR  J*0  TO  Nsets-1 

OUTPUT  @Xy .  Aai  J  , Yal J  > 

NEXT  J 

PRINTER  IS  1 
NE>  T  I 

ASSIGN  9*y  to  * 

SUBEND 
SUE  Njsseit 

COM  /Nus-  I i n , Tsa t , Qdp , Hoc 
Do* .0159 
HO* 1 00C 
IF  I : r»0  --EN 
BEEF 

INcuT  "  ENTER  T5AT  AND  HEAT  FLUX "  ,'sat,  Gdp 
END  IP 

Hf5=PNHf 5' Tsat  ' 

Twc*Tsa t -Qdo /  ho 
Tf=T5at/3»2*Two/3 
K'f-FNK !  T?  ) 

Rho  f  *PNRno( T  f > 

Mu  f =FNMu ( T  f  ) 

Hoc- . BS5*<  Kf '3*Rhof "2»9.81«H*q/< Muf*Do*Qdp  >  >' .333333 
IF  A8S<  <Ho-Hoc  >/Hoc  ).■  .001  THEN 
Ho«<  Ho  +  Hoc  >• . 5 
GOTO  3585 
END  IF 

IF  Im-e  THEN  PRINT  " HO* "  ;  Hoc 

SUBEND 

SUE  Plot  1 


3650 

DIM  v aa >  4 

) 

3655 

PRIN'EP  IS  705 

2EE0 

Idv-1 

3665 

BEEP 

3670 

INPUT  ’OK 

TO  USE 

DEFAULT  VALUES  ( 1 =DEF  =  Y , 0-N  ) " , Idv 

3675 

IF  Idv-1 

THEN 

3680 

Itn-2 

3665 

Xm i n«  1 

3690 

Xnax-5 

3695 

Xstep* 

1 

3700 

Ymin-0 

3705 

Ynax-2 

.0 

3710 

Ystep- 

.5 

3715 

EtSE 

3725 

INPUT 

’ENTER 

MINIMUM  AND  MAXIMUM  X-UALUES" , Xm n , Xma 

3730 

BEEP 

3735 

INPUT 

’ENTER 

MINIMUM  AND  MAXIMUM  Y-UALUES " , Ynm , Yna 

3740 

BEEP 

3745 

INPUT 

"ENTER 

STEP  SIZE  FOR  X-AX I 5 ’ , X s t ep 

3750 

BEEP 

3755 

INPUT 

’ENTER 

STEP  SIZE  FOR  Y-AXIS” , Ys tep 
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END  IF 

PRINT  " IN. SRI  :  IF  230C, 1B00, 552? .5503 . " 


2”-0 

dr  I  NT  "53  0, 100,0, 100' tl 

~  ■y 

3-75 

5  •'  »  «  1 00  ■  '  Xma*-Xmin  .' 

3750 

5 1 y  =  1 03/ f  Yma* - * m 1 n  > 

3"55 

PRINT  "=L  C.C  RD" 

3-90 

FOP  Xa-x- : r  "2  Xma*  STEP 

V  C' 

3~95 

X  » 1  x  a  -  X  rtml'iN 

3500 

PRINT  "  FA C .  XT 

3905 

NE  x  T  x  a 

3510 

PRINT  "PA  100,0. PU.  " 

3315 

PRINT  " P L  PA  0,0  PD" 

3520 

FOR  Ya-Ymin  TO  Yma,"  STEF 

V  s  t  e  p 

»  r E 

Y  = ' Ya-rmin  )•= ey 

3630 

PRINT  "PA  0, " . Y , “YT" 

3835 

NEXT  'i  a 

3640 

POINT  “FA  0,100  Tl.  0  2~ 

3645 

FOR  Xa-Xmir  TO  Xma*  STEF 

x-tep 

3650 

X  =  t  xa-Xm. in/«STx 

3855 

PRINT  "PA'  , X , "  ,  100.  . 

>  7  " 

35E0 

Nf  > T  X  a 

3  555 

PRINT  "PA  100,100  PL  FA 

100.0  DC 

36~0 

f;r  Ya=*m;r  to  >ma*  STEP 

v.;*p 

39-j 

Y= ■  *a-*m:n/*5fy 

3?  50 

FRINT  “P2  PA  10C,  '  ,v 

,  "  Y  T  " 

NEXT  v3 

Pp I NT  ' PA  100,100  F'J' 

—  C  -  'Z 

PRINT  "PA  0,-3  SF  1.9.2" 

35  00 

FOR  Xa-Xmir  TO  Xma*  STEP 

X  s  t  ep 

3905 

X - ( Xa-Xm;n  ;  *  S  Y  * 

3510 

PR  IN’  "PA'  .X,'  ,0.  ' 

3515 

PRINT  *cp  -2,-i .LE" . 

Xa  ;  " " 

3920 

NEXT  xa 

3925 

PPINT  “PU  PA  0,0" 

3930 

-OP  Ye-xmir.  ~0  Yma*  STEP 

Ysteo 

3935 

Y  « .  Y  a  -  y  m :  n  >  *  5  T  y 

3940 

PRINT  -pa  0, ” . Y, "" 

3  945 

PPINT  "CP  -4 , -  .  25  ■  Le 

" , Ya  .  “ " 

3S50  NEx’  V a 


3=5=  BEE- 

3950  I  NP'JT  "SElECT  MODE  (  0«HN/Hi  ,  1»HN<  avg  >/Hl  !",Is"> 

3  3ES  I sm= I sm* 1 

3970  IF  Idv-l  THEN 

397=  IF  Ism-1  THEN  YlabelS- “HN/H1 " 

3990  IF  Ism«2  THEN  Y1 abe IS- ”HN' a vg  )/Hl " 

3995  X 1 abe 1 $=' Tube  Number" 


3990  ELSE 

3995  BEEP 

4000  INPUT  -ENTEP  X-LA8EL" .XlabelS 

4005  BEEP 

4010.  INPUT  "ENTER  Y-LABEL" .YlabelS 

401=  ENO  IF 

4020  PRINT  "SR  1 . 5 , 2 ; PU  PA  50. -10  CP  * . -LEN( X I abe 1 S ) '2 ; "0 ; LB "  ; X I abe  1  * i 
4025  PRINT  "PA  -11,50  CP  0, ' . -LEN< YlabelS >/2*S/B . "DI  0 , 1 . LB' • YlabelS . 
4030  PRINT  "CP  0,0  DI" 

4035  0*0-1 


4040  BEEP 

4045  INPUT  "  WANT  TO  PlOT  DATA  FPOM  A  FILE  (  1  -DEF- Y  ,  0*=N  )'> "  ,  Ok  p 
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4055 

4050 

4055 

4070 

4075 

4000 

4085 

4050 

4055 

4  ]  00 

4105 

4110 

4115 

4120 

4125 

4130 

4135 

4140 

4145 

4150 

4155 

4150 

4165 

4  170 

4175 

4180 

4185 

4150 

4155 

4200 

4205 

4210 

4215 

4220 

4225 

4230 

4235 

4240 

4245 

4250 

4255 

4260 

4265 

4270 

4275 

4280 

4285 

4290 

42  95 

4300 

4305 

4310 

4315 

4320 

4325 

4330 

4335 

4340 

4345 

4350 


BEEF 

INPl|T  "ENTER  THE  NAME  OF  THE  DATA  FI  WE  "  ,  Of  i 1  el 

ASSIGN  @F:le  TO  DfileS 

BEEF 

INPUT  "ENTER  THE  BEGINNING  RUN  NUMBER'  . fid 
BEEF 

INPUT  "ENTER  THE  NUMBER  OF  X-'»  PAIRS  STORED" , Nsets 
BEEP 

input  "select  a  symbol  for  the  plotter  < i=.,2*-,3=c,4= 

FFINT  "PU  DI " 

IF  Sy-1  THEN  PRINT  "SM*" 

IF  Sy-2  THEN  PRINT  "SM*" 

IF  Sy  =  3  THEN  PRINT  "SMc" 

IF  Sv-4  THEN  PRINT  “SMo" 

IF  5  >  =5  THEN  PRINT  "SM"" 

FOR  I  •  1  TO  Nsets 
FOR  J-0  TO  3 

ENTER  @Fi le ; Vaa' J  !,D 
IF  J-0  THEN  Ytop=vaav0! 

Vaat  J  )  -  v  a  a 1  j  i  /  v  t  op 
NE*T  J 

FOR  J-0  TO  3 

>  - !  J  +  1-  Xmj.n  ) ♦  5  f  > 

Y  » ' Y  a  a  i J  i-Ynir  ' • S f y 
PRINT  "PA" ,  X , Y, "" 

NExT  J 
NE*T  I 
BEE  = 

ASSIGN  @Riie  TO  • 

GOTO  4040 
END  IF 
BEER 

INPUT  "LIKE  TO  PLOT  THE  NULSELT  RELATION  <  1-Y , 0-N OKnu 
PRINT  "PU ; SM" 

if  Ox nu«  =  i  then 

FOR  X  a  -  X  m  i  n  TO  Xna*  pT£p  >sted''50 
X -t  Xa-Xmi n  ) • S  f a 

IF  Iam-I  ANO  Xa'-Xrun  THEN  Ya-Xa "  . 75- ■;  X a- 1  : "  . 75 
IF  Ian-2  AND  Xa  >Xmjn  then  Ya-Xa" t - . 25  ) 

IF  Xa-Xmin  THEN  Ya-1 
Y-t Ya-Ynin >«Sfy 
PRINT  "PA' , X , Y , *PD" 

NEXT  Xa 
BEEF 

PRINT  "PU" 

INPUT  "MOUE  THE  PEN  TO  LABEl  THE  NUSSEuT  LINE", OK 
PRINT  "LBNussel t  " 

END  IF 

IF  Ism-2  then 
BEEP 

INPUT  "LIKE  TO  PLOT  EXPTl  CURUE  (  1  -Y  ,  0-N  ) "  ,  OKex 
Nq-0 

IF  OK  e>  *  1  THEN 
BEEP 

INPUT  "ENTER  THE  EXPONENT", Ex 
FOR  Xa-Xmin  TO  Xmax  STEP  Xstep/10 
Nq-Nq+1 
Ya-Xa'(-Ex  ) 

X»( Xa-Xmin  >*Sf  x 
Y»< Ya-Ymin  >«Sfy 


36 


435= 
4353 
435  = 


43-  = 


438  = 


IF  N3  MOD  Z-0  THEN 

CR I  NT  “ PA " , X , Y  ,  CD" 
ELSE 

OR  I  NT  1  PA "  ,  *  .  Y  , “CU ' 
END  IF 
ne  *  >  a 

PRINT  'On" 


435=  INF.-  "  MQL’E  PEN  TO  lASEl  AND  HIT  ENTER "  ,  3  x 

4*20  PRINT  " L E s “ 2 " 

4425  PRINT  "PR  -1  0" 

44  12  PRINT  * L E " . E * • "  " 

4415  GOTO  4302 

4432  END  IF 

4435  END  IF 

4432  GOT'  4512 

4435  EEE= 

4442  INPlt  "  1_  I H  E  T3  Plot  fern  RELATIONSHIP  i-Y  ,2=N  n- .  Yes 
4445  IF  Yes-1  *HEN 

4452  F  OF  x a* * n x p  t n  xnqa,  J'EP  X  stec  -’30 

4455  >a-«a' f-1  6  > 

4452  x  =  <  *  a  - Xp :  n  1  *  5  *  / 

445S  «'< ia-'pin  1  ♦ S  f  > 

4472  POINT  'FA'  , x , o ,  "or • 

44~5  NEXT  Xa 

4482  PR  IN’  "  =  U" 

4465  BEE- 

4452  INPlt  "MOLE  THE  PEN  TC  LABEL  KERN  RElAT I ONShI F ' , 0> 

4495  PRINT  " LBKepp : P J " 

4S22  END  IF 

4525  PRINT  "PU  PA  2,0" 

4512  EEE= 

451=  INRU*  "lIKE  TC  PlO”  EISSENBER6  RELATION  ( i-y,0»N  , OX ei 

4532  IF  Cx  e  1  =  I  THEN 

453=  FOR  xa  =  x,(r;n  TO  Xna*  £TpF  Xstep/10 

4532  xa=  .£♦ .43*Xa '  -.35) 

4  535  x«(Xa-Xi*>io  •  5 f > 

4543  y« <  Ya- ’>.  1  p  > •  S f  > 

454=  PRINT  PA" ,X  ,X  , "CP- 

4552  NEXT  X a 


4555  PRINT  "FU" 

4552  EEE= 

45E5  INPLT  " MOl'E  THE  PEN  TO  LABEL  THE  EIS5ENBERG  LINE", Ok 

4572  PPINr  " LBE : s sencepg ■ Pj " 

4575  END  1= 

4582  PR IN7  -p-j  s P0" 

4585  SUBENu 

4590  DEF  FNPvst ( Tc  ) 

4SSS  COM  /Fid/  Ift 

4B02  DIM  K C  8  ) 

4605  IF  Ift-0  THEN 

4610  BEEF 

461=  PRINT  "POST  CORRELATION  NOT  AVAILABLE  FOR  R-114” 

4E20  STOP 

4625  END  IF 

4630  IF  Ift-i  THEN 

4635  DATA  -7.691334564,-26.08033696.-168.1706646,64.23285504,-116.9646 

4640  DATA  4. 16711732,20.9750675,1 .£9,6 

4645  PEAD  K< •  > 
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4650  T-(Tc  +  273. 15  1/647.3 

4655  Sum-0 

4660  FOR  N-0  TO  4 

4665  Sum-Sum+K(N  )•( 1-T  )‘<N+1  ) 

4670  NEXT  N 

4675  Br-Sum/(  T*<  1+K  ( 5  >*  (  1-T  )+K(  6  )•  (  1-T  )*2  )  >-(  1-T  )/{  K  ( 7  )•  (  1-T  )*2+K(  8  )  ) 

4680  Pr-EXP(Br) 

4665  P-22l20000*Pr 

4690  END  IF 

4695  IF  Ift-2  THEN 

4700  Tf-Tc»l .8+32+459.6 

4705  P-10'(33. 0655-4330. 98/T f-9 . 2635»LGT< T f  )+2 . 0539E-3+T f  ) 

4715  END  IF 

4720  IF  Ift-3  THEN 

4  725  A-9. 394685-3066. l/( Tc  +  273. 15  ) 

4730  P-133. 32*10‘A 

4735  END  IF 

4740  RETURN  P 

4745  FNEND 

4750  DEF  FNHf g( T  ) 

4755  COM  /Fid/  Ift 

4760  IF  Ift-0  THEN 

4765  T  f »T • 1 . 8+32 

4770  Hfg-6. 14S15S8E+1 -Tf ♦<  6 . 951079E-2+T f • ( 1 . 3988686E-4+ 1 . 9607843E-7. t f 

4775  Hfg=Hf g+2326 

4780  END  IF 

4785  IF  Ift-l  THEN 

4790  Hfg-2477200-2450.(  T-10 ) 

4795  END  IF 

4800  IF  Ift  —2  THEN 

4805  T  f-T»  1 . 8  +  32 

4810  Hf g»7 . 05S7857E+1-T  f •( 4 . 838052E-2  +  1 ,2619048E-4»Tf > 

4815  Hfg-Hf g*2326 

4820  END  IF 

4825  IF  Ift-3  THEN 

4830  Tk-T+273.15 

4835  Hfg-1 . 35264E+6-Tk»( 6 . 38263E+2+TM . 747462 > 

4840  END  IF 

4645  RETURN  Hf g 

4850  FNEND 

4855  DEF  FNMu( T  ) 

4850  COM  /Fid/  Ift 

4855  IF  Ift-0  THEN 

4870  Tk-T+273.15 

4875  Mu-EXP(-4. 4636+1011 ,47/Tk  l+l.E-3 

4880  END  IF 

4885  IF  Ift-l  THEN 

4890  A-247.8/1 T+133. 15 ) 

4895  Mu«2.4E-5»10'A 

4900  END  IF 

4905  IF  Ift-2  THEN 

4910  Mu-8. 962981 9E-4-T+< 1 . 1094609E-5-T+5.566829E-8 ) 

4915  END  IF 

4920  IF  Ift-3  THEN 

4925  Tk- 1 / < T  +  273. 15  ) 

4930  Mu-EXP ( -11 . 0179+Tk  *  <  1 . 744E  +  3-Tk • ( 2 . 80335E  +  5-Tk • 1 . 1266 1E  +  8 > ) ) 

4935  END  IF 

4940  RETURN  Mu 

4945  FNEND 
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I 


4950 

4955 

4960 

4965 

4970 

4975 

4980 

4985 

4990 

4995 

5000 

5005 

5010 

5015 

5020 

5025 

5030 

5035 

bk)4tf 

5045 

5050 

5055 

5060 

5065 

5070 

5075 

5080 

5085 

5090 

5095 

5100 

5105 

5110 

5115 

5120 

5125 

5130 

5135 

5140 

5145 

5150 

5155 

5160 

5165 

5170 

5175 

5180 

5185 

5190 

5195 

5200 

5205 

5210 

5215 

5220 

5225 

5230 

5235 

5240 


DEF  FNVvst ( T t  ) 

COM  /Fid/  1ft 
IF  lft«0  THEN 
BEEP 

PRINT  -UO ST  CORRELATION  NOT  AVAILABLE  FOR  R-114* 

STOP 
END  IF 

IF  Ift-1  THEN 
P-FNPvstl Tt ) 

T-Tt+273. 15 
X-1S00/T 

Fl-l/i 1+T+ 1 . E-4  ) 

F2-< 1 —EXP  (  -X ) )'2.5+EXP(X  )/X*.5 
B-.001S+F1-. 000942 +F2 - . 0004882 »X 
K-2»P/<  461 .52  +  T  ) 

V-<  l  +  (  1+2+B  +  K  )*  .5  >/K 
END  IF 

IF  Ift-2  THEN 
I t- 1 1  *  1 . S  +  Jt 

V*13 . 955357-T  f  • ' . 161272E2-Tf*S. 1726190E-4 ) 

V-V/16.018 
END  IF 

IF  Ift-3  THEN 
Tk-Tt+273.1S 
P-FNPvst i  Tt  ) 

V-133 . 95+Tk/P 
END  IF 
RETURN  V 
FNENO 

DEF  FNCp(T) 

COM  /Fid/  Ift 
IF  Ift-0  THEN 
Tk-T+273.15 

Cp” . 401 1 8+Tk • ( 1 .65007E-3+Tk+( 1 .5l4  94E-6-Tk»6 .67853E-10  >  ) 

END  IF 

IF  Ift-1  THEN 

Cp-4 . 2 1 1 20858  — T * ( 2.26826E-3-T»<4.4236lE-5+2  71428E-7+T)! 

END  IF 

IF  Ift-2  THEN 

Cp-9.2507273E-1+T*( 9 . 3400433E-4 + l . 7207792E-6 »T > 

END  IF 

IF  Ift-3  THEN 
Tk-T+273.15 

Cp-4. 1868*  (  1.6884E-2  +  Tk»(3. 35083E-3-TI.  •  (  7 . 224E-6-Tk  *7 . El 748E-9  )  > 
END  IF 

RETURN  Cp+1000 
FNEND 

DEF  FNRho< T  ) 

COM  /Fid/  Ift 
IF  Ift-0  THEN 
Tk-T+273.15 
X - 1 - ! 1 ,8»Tk/753.95  ) 

Ro-36. 32+61 .146414  +  X"(  1/3  I+16.418015+X+17.476838+X- .5+1 .  1 1 ge^B* x ' 
Ro-Ro/ . 062428 
END  IF 

IF  Ift-1  THEN 

Ro-999 . 52946+T+ ( . 01 269-T* ( 5 , 48251 3E-3-T+ 1 .234147E-5 ) ) 

END  IF 

IF  Ift-2  THEN 
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5245  Rc-1 .eZe“4-79E-Z-T*<  2.Zi8B346-T.2.2ET5251E-3 > 

52=0  END  IF 

52EE  1  =  I  f  *  =3  THEN 

52E0  Tk = '*273 . 1 5-326 . 1 5 

5255  Uf-9.2484SE-4+Tk«<5. 2''96E-7"T>  •iS.2444E-10+T**3.057E- 

5270  Ro= 1 /V f 

52“5  END  IF 

52 SC  RETURN  Ro 

5285  FNEND 

5290  DEF  FNPr( T  ) 

5295  Fr-FNCp' T  i  +  FNMut T '/FNK< T > 

5300  RETURN  Pr 

S3C5  FNEND 

5310  DEF  FNK  f  T  i 

5315  COM  /Fid/  Ift 

5320  IF  Ift-0  THEN  K- . 071 00026 1  * T 

5325  IF  Ift-1  THEN 

5330  y-' T  +  2'73 . 15  V'2',3 . 15 

5235  K-- . 92247+X • < 2 .  839S-X»<  1 ,8007-X*<  .  52577- . 07344  +  X  )  )  ) 

5340  END  IF 

5345  IF  I f t -2  THEN 

5350  K“6.2095238E-2-T»( 2 . 22 1428SE-4+T+2 . 2809524E-8 > 

5255  END  IF 

53EC  IF  I  f t-3  THEN 

53E5  Tk-T+273.15 

52'5  END  IF 

5380  RETURN  k 

53E5  FNEND 

5350  DEF  FNHf(T) 

5395  COM  /Fid/  Ift 

5400  IF  Ift-0  THEN 

5405  BEEP 

5410  PRINT  "HF  CORRELATION  NOT  FOR  R-114" 

54.5  STOP 

5420  END  IF 

5425  If  I't-i  THEN 

5430  Hf-T+<4.203849-T«'  5.88122E-4-T.4.55160317E-6  >: 

5435  END  IF 

-.440  IF  I  f  t * 2  Then 

544=  T  f -t. i . g»32 

-452  h  f  »S  .  GCSST 1  +T  f  •  .194E”85'  +  Tf»i.3214285E-4; 

5455  hf -Hf »2 . 32E 

-4=0  EN2  I f 

54=5  IF  ift-3  Then 

=4-0  Ht«250  ■  TO  BE  VERIFIED 

54-5  END  IF 

5480  RETURN  Hf  *  1000 

5485  FNEND 

5450  SUB  Plot2 

5495  COM  /Drl/  S t ar , Sym , I  con 

5500  COM  /Fid/  Ift 

5505  DIM  Ci9),Xya<7),Doa(3) 

5510  DATA  0.0158,0.0158,0.0156,0.0158 

5515  RE AO  Doa< •  ) 

5520  Fu-1 

5525  PR  INTER  IS  1 

5520  8EEP 

5535  PRINT  USING  ”  4X  ,  “ "  Se  1  ec  t  Ootion  X-Y  Units 

5540  PRINT  USING  "Ex,""0  Use  de'ault  values 
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554  = 
5550 
5555 
5560 

55-0 

5=-S 

5580 

c  =  =  = 

SS9C 

5595 

560C 

5605 

5610 

55:5 

5622 

5525 

5630 

5535 

5640 

5545 

5650 

5655 

5660 

5665 

5670 

5575 

5680 

5635 

5690 

5655 

5  .'00 

570= 

5710 

5715 

5~2C 

5725 

5720 

5’25 

5740 

5~45 

5”5C 

5755 

5750 

5-65 

5"70 

S~7S 

5780 

5785 

5790 

5795 

5800 

5805 

5810 

5815 

5820 

5825 

5820 

5855 

5840 


c P 1  n T  _c;mc-  use  net  values 

I 2- d 

BEEF 

INPUT  “  E  N  T  E  c  TUBE  C3DE " , I  cede 
Do“Dca -  I  tube 
Iht  =  2 
BEEP 

PP1NT  USING  "4* , “"Select  option  """ 

PRINT  USING  “6>.  .'"'0  n  ven«us  c  ' 

PRINT  USING  “  E  X  , "  “  1  q  versus  Delta-T""' 

PRINT  USING  "6X,'"2  h  versus  Oelta-T  (default) . 

INPUT  I bt 
PRINTER  IS  705 

IF  Okd-0  THEN  (AXIS  DEFAULT  VALUES 
IF  Iht-0  THEN  '  <  h  vs  g i 
Ymn»0 
Yma’,*60 
Ystep-10 
Xnin* . 2 
Xnax* 1 . 4 
Xstep*.2 
END  IF 

Ir  lHt-1  THEN  1  (  g  vs  t  ) 

Xmin»0 
Ynin«0 
Ymax« . 5 
Xmax«15 
Xstep-3 
Ys  tep  = . 1 
END  IF 

IF  Int-2  THEN  ' ( h  vs  t  ) 

Yni  n»k5 
Yni n*0 
Xna* ”SC 
vna>  *6 
Xstep=10 
v  s  t  ep  =  1 
END  IF 

END 

IF  Of.  d*l  THEN 
BEEF 

INPUT  '  EN'EF  MINI  HUM  AND  MAXIMUM  X  -UAlUES  "  ,  Xn.ir  ,  »(^a» 
BEEF 

INPUT  “ENTER  MINIMUM  AND  MAXIMUM  Y-VAL'JES" , Yrin, Yma* 
BEEF 

INFU"  “ENTER  STEF  SIZE  FOR  X -AX 1 5 " , Xs t ep 
BEEF 

IN=>Ut  “ENTER  STEP  SI2E  FOR  Y-ftx IS"  ,  Ys t ep 
END  IF 
BEEF 

PRINT  " I N ; 5 P 1 . I P  2300, 1800,8300,6800. " 

PP I NT  “SC  0,100,0, 100. TL  2,0;" 

Sf»«100  "  Xna.r-Xnin  ) 

S  f  y* 1 00  / ' Ynax-Ynin ) 

BEEP 

Icq«0 

INPUT  ‘LIKE  TO  BY-PASS  CAGE  ( 1  - Y , 0-N-DEFAULT P ‘ , I  eg 
IF  Icg*l  THEN  6175 
°P  I  NT  'PU  0,0  FD“ 

FQF  XaMXm  1  r.  TC  Xna  -  STEF  Xstep 
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5845 
5853 
5855 
5880 
5365 
5570 
537  = 
5330 
5885 
5890 
5895 
5900 
5505 
5910 
5915 
5920 
5525 
5930 
5935 
5940 
5945 
5950 
5955 
5980 
5965 
5970 
5975 
5980 
5985 
5990 
5995 
6000 
6005 
6010 
6015 
6020 
6025 

6035 
6040 
604  = 
6050 
6055 
6060  • 
6065 
6070 
6075 
6080 
6085 
6090 
6095 
6100 
6105 
kUI/n" 
6110 
6115 
6120 
R  1.0, 
6125 
6130 


X  =  (  X. a-Xnm  i » S  f  > 

CRIN*  "PA" ; v , ■  , 0 .  XT;  * 

NEXT  Xa 

PRINT  “FA  100.0.F..‘.  - 
PRINT  ”P'J  PA  0.0  PD" 

FOP  'a=Ymin  TO  fna.-.  STEr  Ystep 
Y» ;  Ye-Ycn  n  ) *S  f  y 
PRINT  "PA  0, ",Y. "VT" 

NEXT  Ya 

PRINT  "PA  0,100  TL  02" 

FOR  Xa*Xrun  TO  Xma*  STEF  Xstep 
X*(  Xa-Xmin  (•Sf* 

PRINT  "PA" ; X , " , 100,  XT" 

NEXT  Xa 

PRINT  ‘PA  100,100  PU  PA  100,0  PO" 

FOP  Ya=Yrun  TO  Yma*  STEF  Ystep 
Y  =  (Ya-Yrun  )»Sfy 
PRINT  "PD  PA  100,  "  ,Y  ,  "YT" 

NEXT  Ya 

PRINT  "Oft  ]00,I00  PU" 

PRINT  "PA  0,-2  SR  1.5,2" 

FOR  Xa-Xmin  TO  Xma*  STEP  Xstep 
X-( Xa-Xmin ) • 5  f  x 
PRINT  " PA " ; X , " ,0,  " 

IF  Xa<  1  AND  Xa-.  >0  THEN  PRINT  "CP  -1 . 5  ,  - 1  .  LB0 ;  PR  -1 , 0 ;  LB  “  .  Xa  ;  "  " 

IF  Xa-0  THEN  PRINT  "CP  -.5.-1.LB0" 

X  i  n-0 

IF  Xa  MOD  1-0  THEN  Xin-1 

IF  Xa>-10  THEN  PRINT  "CP  -2 , - 1 ; LB" . Xa . ’ * 

IF  Xa>l  AND  Xa< 10  AND  Xin-1  THEN  PRINT  "CP  - 1 . 25, -1 ; LB' ;  Xa ;  ” " 

IF  Xa>l  AND  Xin-0  THEN  PRINT  "CP  -2  , - 1 , LB " ; Xa , * " 

IF  Xa-1  THEN  PRINT  "CP  -1,~1;LB1.0“ 

NEXT  Xa 

PRINT  "PU  PA  0,0" 

FOR  Ya«Ymin  TO  Ymax  STEP  Ystep 
Y«(Ya-Ymn)oSfy 
PRINT  "PA  0,  "  ,  Y,”" 

IF  Ya-.l  AND  Ya<>0  THEN  PRINT  "CP  -4  ,  -  .  25  .  lB0  ;  PR  -2 , 0 .  LB "  .  Ya  ;  " " 

IP  Ya-0  THEN  PRINT  "CP  -2.-.2S.LB0" 

IF  Ya/1  AND  Iht<2  THEN  PRINT  "CP  -5 . - . 25 . LB " ; Ya . " “ 

IF  Ya/9  AND  Iht-2  THEN  PRINT  "CP  -4 , - . 25 ; UE “ . Ye . “ " 

IF  Ya>0  AND  Ya< 10  AND  Iht-2  THEN  PRINT  "CP  -3 , - . 25 . LB " . Ya ; " " 

IP  Ya-1  THEN  PRINT  "CP  -4 , - . 25 . LB  1 . 0 “ 

NEXT  Ya 
IF  Okd-1  THEN 
BEEP 

INPUT  "ENTER  X-LABEL" .Xlabel* 

BEEP 

INPUT  "ENTER  Y-LABEL” .Ylabell 
END  IF 

IF  I h t < : 1  THEN 

PRINT  "SR  1 . 5 , 2 ; PU  PA  -12.3S.DI  0,l.LBh:PR  1.0.S.LBo;PR  -l,0.5;LB/( 

PRINT  "PR  -1,C.S;SR  1,1.S;LB2.SR  1.5, 2. PR  .5, .5, LB., PR  ,5,0;LBK>" 
ELSE 

PRINT  "PA  -12,39.01  0,1; LBq /<  MW/n ; SR  1,1.5;PR  - 1  , 0 . 5 . LB2 ; SR  1.5.2.P 
0 . 5 , LB  ) " 

END  IP 

IF  Iht-0  THEN 
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E135  PRINT  -SR  1.S.2.PU  PA  40.-10.QI ,LBq/(MW/m;PR  0.5.1.SP  1 . 1 .5, lb: . s= 

1.5,2. PR  .  5 ,  - 1 .  L  B  ) " 

6140  Ei.SE 

5145  ORINT  "DI  PA  3S,-:0,LB(T,PP  ,5,-l.LBs.PR  .S.'.LB-T.PP  -2.4,3  PD  c= 

2,0  Pu.PR  .5,-4,'' 

6150  PRINT  "LBuo.PR  .S,1;LB>/K” 

E 1 55  END  IP 

E 160  PRINT  'CP  0.0  DI" 

E1E5  Xlg«l.£+S 

6170  Xug*- 1 . £+6 

5  i  "75  Xal-50 

S1S0  Yal-95 

6165  Nrun-0 

6190  BEEP 

6195  INPUT  "WANT  TO  PLOT  DATA  FROM  A  FILE  (1 -Y , 0-N  )f " , Ok 

6200  Xll-l.E+6 

6205  Xu  1 •- 1 . E+6 

6210  Okp-0 

6215  IF  Ok-1  THEN 

£220  BEEP 

S225  INPUT  “ENTER  THE  NAME  OF  THE  PlOT  DATA  FILE" , D_f l 1 e$ 

E230  ASSIGN  8Fi  le  TO  D_fi’ei 

E235  IF  Icortb<>0  THEN  6265 

6240  Sx-0 

6245  Sy-0 

5250  Sx2-0 

6255  Sx  y *0 

6260  Md-1 

6265  BEEP 

6270  INPUT  “ENTER  THE  BEGINNING  RUN  NUMBER  (DEF-l>",Md 

6275  Npaira-9 

6280  BEEP 

6285  INPUT  "ENTER  THE  NUM8ER  OF  X-Y  PAIRS  STORED  <DEF-9>“ .Npairs 

6290  Nrun-Nrun+Npairs 

6295  PRINTER  IS  1 

6300  BEEP 

E305  PRINT  USING  "4X, •“Select  a  symbol  "““ 

6310  PRINT  USING  “4X,““  1  Star  2  Plus  sign""" 

6315  PRINT  USING  "4X,-“  3  Circle  4  Square""" 

6320  PRINT  USING  "4X,"“  5  Rombus“““ 

6325  PRINT  USING  “4X,”"  6  Right-side-up  triangle'"" 

6330  PRINT  USING  "4X,““  7  Up-s i de-douin  triangle""" 

6335  INPUT  Sym 

Bfctr 

6345  INPUT  “ENTER  TUBE  NUMBER  FOR  PLOTTING  ( 0-TOP , 1 -SECOND _ )",Itube 

6350  PRINTER  IS  705 

6355  IF  Syn-1  THEN  PRINT  “SM+" 

6360  IP  Syn-2  THEN  PRINT  "SM+" 

6365  IF  Syn-3  THEN  PRINT  “SMo" 

6370  IF  f1d>  1  THEN 

6375  FOR  1-1  TO  ( Md-1  ) 

6330  ENTER  9File;Xya<»  ) 

6395  NEXT  I 

6390  END  IF 

6395  FOR  1-1  TO  Npairs 

6400  ENTER  0Fi le ; Xya<  «  ) 

6405  Ya-Xya( I tube»2  ) 

6410  Xa-Xya( I  tube-2  +  1  ) 

6415  Yc-L0G( Xa  ) 
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uj  uj  in 


£425 

430 

440 
644  = 

6450 
6455 
E46C 
6465 
£470 
6475 
£480 
6485 
6450 
6495 
£500 
650  = 
£510 


£50  = 
£530 
6535 
£540 
£545 
6550 
6555 
6560 
65ES 
5570 
6575 
6560 
6585 
£550 
6565 
6600 
6S0S 
6510 
6515 
£620 
6625 
6500 
6625 
6E40 
6645 
6650 
6655 
6660 
6665 
6675 
6680 
6685 
6690 
66  95 
6700 
6705 
6710 
£715 


Xc-lOS; >6' >a  ) 
S *  -  5*  xx  c 
Sy-Sy+Yc 
S*2«5>  2  +  >  c ' 2 

£,y«S*>-r.«C*>C 

IF  Iht-0  THEN 
Xt-Ya 
Ya*ra/Xa 


Xa  =  xt 

IF  Xa/1 .E*6>Xjl  THEN  Xui-Xa/l.E+6 
IF  XeM.E+B'Xli  THEN  Xll-Xa/l.E+6 
END  IF 

IF  Iht *0  THEN 

X»<  Xa»l  .E-E-Xmin  ) •  S f  >■ 
y  «  i  y  a  *  1 .  E-  3-  Yn  r  >  •  S  f  y 
END  IF 

IF  Irt=l  THEM 

y«i  x a-Xmn  )»S f ' 

Y«<  Va«l  .E-E-Ymr.  ;  *Sf  y 
END  I? 

If  I  r,  t  =  I  THEN 

X  *  ■  >  a  /  Y  a  -  xr :  n  f  • 

y«< Ya*l  .E-3-'n:r  ;*5‘y 


END  IF 

IF  Y >  1 00  OR  Y-:0  THEN  E555 
IF  Syf»-'3  THEN  PRINT  "SM" 

IF  Syn<4  THEN  PRINT  "SR  1.4, 2. 4" 

PRINT  ”  P  A "  ,  X  ,  Y  ,  " 

IF  Sym>3  THEN  PRINT  "SR  1.2, 1.6" 

IF  Syn-4  THEN  PRINT  " UC2 , 4 , 99 , 0 , -8 . -4 , 0 , 0 . 8 . 4 . 0 ; ” 

IF  Syn-5  THEN  PRINT  " UC3 . 0 , 99 , -3 , -6 , -3 , 6 . 3 , 6 . 2 , -6 ; * 
IF  Sym-6  THEN  PRINT  “ UC0.S . 3 . 99 , 3 , -8 . -6 . 0 , 3 . 8 . " 

IF  Syn-7  THEN  PRINT  "UC0, -5. 3, 99, -3, 8, 6,0, -3, ~8;" 


NEXT  I 
BEEP 

INPUT  "WANT  TO  LABEL  <  1-Y ,0*N  ,  I  lb! 

IF  Ilbl-1  THEN 

IF  Syn>3  THEN  PRINT  "SM" 

1=  5yH<  4  THEN  PRINT  "SR  1  .  4 , 2 . 4 " 

PRINT  "PA" ,Xal . Yal , 

IF  Synx3  THEN  PRINT  "SP  1.2, 1.6" 

IF  Syn-4  THEN  PRINT  "UC2 , 4 , 99 , 0 , -8 , -4 , 0 , 0 , 8 , 4 , 0 ; ' 
IF  Syn-5  THEN  PRINT  "UC3 , 0 , 99 , -3 . -6 , -3 , 6 , 3 , 6 , 2 , -6 ; 
IF  Syn-6  THEN  PRINT  "UC0 , 5 . 3 , 99 , 3 , -8 , -6 , 0 , 3 , 8 ; “ 

IF  Syn-7  THEN  PRINT  'UC0 , -5 . 3 , 99 , -3 , 8 . 6 , 0 , -3 , -8 , ‘ 
PRINT  "SM" 

IF  Syn<4  THEN  PRINT  "PR  2,0" 

PRINT  "PR  2 ,-l .0; SR  1 . 0 , 1 . 8 ; LB" ; D_f 1 le* . " ” 
Yal-Yai-S 


BEEP 

IF  Ias-1  THEN 
BEEF 

INPUT  "ENTER  THE  STRING" .Label* 

PRINT  'PR  2,0, SR  1. 0,1. 8. LB". Label*;"" 
GOTO  6665 
END  IF 
END  IF 
BEEP 


INPUT  "WANT  TO  COMBINE  ANOTHER  FILE’  < 1 -Y , 0-N  ) " , I conb 
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6720 

6725 

6730 

6735 

6740 

6745 

6750 

6755 

6760 

6765 

6770 

6775 

6780 

6785 

6790 

6795 

6800 

6906 

6910 

6615 

6620 

6525 

6530 

6635 

6S4C 

6545 

6850 

6655 

6860 

6865 

6870 

6875 

6880 

6885 

6830 

6895 

6900 

5,05 

6910 

6915 

6920 

6925 

6930 

6935 

6940 

6949 

6950 
6955 
6960 
6965 
6970 
6975 
6980 
6985 
6990 
6995 
/MV 
7005 
7010 


AS5I6N  fEjl4  TO  » 

X12-5 

Xul«4S 

IF  IconO <  .>0  THEN  6220 
BEEP 

input  -want  to  plot  a  least-squares  line  < i-y.o-nt . iu 

IF  Ils«l  THEN 
BEEP 

INPUT  ‘ SELECT  EXPONENT :  0-COMPUTE,  1 -0 . 75 > " , lexp 
BEEP 

INPUT  'SELECT  CURU6  TYPE  (0-SOlID. 1-DASHED  )*. I  It 

IIt-Ilt+1 

PRINT  -SM- 

IF  Iexp-0  THEN 

8b-(Nrun*Sxy_Sy*Sx  )/(Nrun*Sx2-Sx"2  > 

ELSE 

6b-. 75 
EN0  IF 

Aa=< Sy"Bb*Sx  )/Nrun 
Aa-EXPi Aa  ) 

PRINTER  IS  1 

PRINT  USING  “lex.-'a  -  "",Z.4DE".Aa 
PRINT  USING  " 1 0 X , " " n  -  '",Z.4D£"  6d 
PRINTER  IS  705 
In-0 

IF  Iht-0  THEN  Xx»lep*xstep.'40 
IF  lht>0  THEN  Xxstep-Xstep/ 10 
FOR  Xa-XU  TO  Xu  I  STEF  Xxstep 
IF  Xa>.99«XKi0x  THEN  6995 
IF  IHt-1  THEN  Ya-Aa»Xa “Bb 

IF  iht-0  THEN  Ya-Aa'<  1/Bb  >•(  Xa*l  .E*6  >'(  (Bb-n/Bb  ) 
IF  I  a  t  *2  THEN  Ya-Aa*Xa“<  Bb-t  ) 

IF  Iht-0  THEN 

Y»(Ya*l.E-3-Ymin)«Sfy 
X“(  Xa-Xrnin  )*Sf  x 
END  IF 

IF  Iht-1  THEN 

Y-<Y««l.E_6-Yr>iin)»Sfy 
X-( Xa-Xhin  )*Sf x 
END  IF 

IF  Iht-2  THEN 

Y-(Ya*l.E-3-Yi'un)*Sfy 
X-( Xa-Xnin  )»Sf x 
END  IF 

IF  Y <0  THEN  Y-0 
IF  Y> 100  THEN  6990 
IF  Ilt-1  THEN 

PRINT  'PA' ,X,Y, -PO" 

ELSE 

In«In+ 1 
Ir-In  MOD  Ill 

IF  Ir-1  THEN  PRINT  "PA* ,X , Y , *PD" 

IF  lr-Q  THEN  PRINT  "PA" , X , Y, "PU* 

END  IF 
NEXT  Xa 
PRINT  *PU' 
tNU  ih 
Iconb-0 
60T 0  SI8S 
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END  IF 

PRINT  *P'J  SM " 

SEEP 

input  -want  to  plot  nusselT  line  <  i-»  ,e-N  . 

IF  lnp-0  Then  7125 

8EEP 

INPUT  "ENTER  TS  AT  (DEFAULT-18  DEG  C!",Tsat 
h'g-FNHFg'.  T sat  ) 

>11-5 

Xul*45 

FOR  Xa-Xll  TO  X u I  STEP  Xstep/SO 
T  f 1 1 m-Tsa t -X  a* .5 
K  f-FNK ( T  f l In > 

Rhof -FNRhoi  T  f  i  If".  / 

Muf-PNMu<  Tf iln ) 

Ya-.728*(Kf"3,RHof"2#9.81«H>g',<Muf»Do»>a 
X-<  Xa-Xmin  )»Sf  » 

Y - ( Y  a  •  1 .E-3-Ynjr i  •  S  f  y 
PRINT  "PA" ,X , Y. "FD” 
nE  x  T  x  a 

PRINT  "PU  PA  0,0" 

PRINT  "P'J  PA  0,0  SP0 " 

SuBENO 
SUE  p;3t 3 

COd  /Dr-1/  S  t  ar  ,  S  ,  n  .  Iccr, 


"140 

007  /Fid/  Ift 

7145 

cm  ct  9 . 

"150 

Fu*  1 

7155 

PRINTER  IS  1 

-;S0 

BEEP 

7155 

PRINT  USING  *'4X 

"170 

PRINT  USING  "GX 

7;  7= 

PR  IN"  USING  "Gx 

T  1  S0 

INPut  o> a 

7185 

PRINTER  IS  70S 

"1  S0 

if  o> a-0  then 

71S5 

Xmi n-0 

"200 

Ynir*0 

"205 

Xna»-  1 5 

"210 

Yna-*15 

7215 

Xstep-3 

7200 

Ystep-3 

"Select  Opticn  X-Y  Units 
"0  Use  default  values'"" 

"  I  Use  new  values"' 


“’230 

BEEP 

7235 

INPUT 

"ENTER  MINI-JM  AND  MAXIMUM  X-UAlUES'.Xmn.Xnax 

7240 

BEEP 

7245 

INPUT 

"ENTER  MINIMUM  AND  MAXIMUM  Y-UALUES " , Ym n , Ynax 

7250 

BEEP 

7255 

INPUT 

"ENTER  STEP  SIZE  FOR  X-AX I S ' . Xs t ep 

7260 

BEEP 

7265 

INPUT 

"ENTER  STEP  SIZE  FOR  Y-AX IS ” . Ys t ep 

7270 

END  IP 

7275 

BEEP 

7280 

PRINT  "IN 

; SP1 ; IP  2300,1800,8300,6800;" 

7285 

PRINT  'SC 

0, 100,0, 100. TL  2,0: ' 

7290 

Sfx-100/< Xmax-Xnin ) 

7295 

Sfy*100/(Ynax-Ymin> 

7300 

BEEP 

7305 

Icg-0 
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'  -  %. 

INPUT  "lIkE  TO  By-PASS  CAGE  (l-Y 

,0*N*O£ FAulT  '/ 7 " ,  Icq 

~  ■  - 

IP  I c o “ 1  Then  7620 

_z:e 

PRINT  'P'J  0,0  PD" 

~  2  2^ 

r3P  x  a*  xh. ; n  TC  xna  ■  S'EP  »  s  t  ec 

==:n*  •=*■.>,' .e  >- . 

“34? 

NEx  '  >  i 

“345 

-PIN*  ■  =  A  100,0.  PI¬ 

“35? 

PPIN"  -pij  pa  0,0  pc 

-?p  >a*v«-ir.  TC  STEP  'stec 

“3E? 

Y  *  >  > a- » m  1  r  ’•ST, 

"255 

pp  IN"  "PA  "  vT" 

NE  >T  > a 

-  7  — - 

=  =I\"  " p A  0,10?  2  c - 

“35? 

-  -  -  *  a*  »  »  :  r  T  D  x  r*t  a  *  STEP  x  ,  ^ 

“35  = 

»  ■ ;  x  a-  •  >• :  n  •  E  £> 

CP I “ P A '  . x , '  , ;0? .  *  - • 

“3  9S 

NE  x  *  t a 

"4?? 

p  P  I  N"  P  a  100,10?  p  ph  i??.?  c?' 

“4?  = 

p C P  ’a*x"*:r  TC  Y>a  >  S"EP  yite^ 

~  4  •  *» 

**  •  V. 

X  »  '  •fi-’ff'ir  •  S  f , 

“4 :  5 

P=IN"  ' P"  PA  ;??, 

4  2  2 

NE'"  •* 

“435 

==:*,-  'pa  >??, loo  pj’ 

~  4  7  3 

PR  IN”  PA  ?- ;  ;c  x.p,; 

“425 

•  —  P  Xaa  x*  i  r  _  X rxfl,  S  pP  1  ,  i gr 

4  4  •» 

**  «•  t 

>  ■ '  '  a  -  *  -  ;  -  ■  •  S  f  « 

744  = 

PR  IN'  "PA' . > , " ,0. " 

“45? 

IP  xa-  1  AND  x  a  ■  >0  THEN  =PIN" 

•  ;p  -1.5.-;. ug?  pc  - ; 

“455 

;p  xa  =  ?  Then  PRINT  "CP  -.5,-1 

...60" 

“45C 

X  i  n«0 

“455 

IP  xa  MCC  1*0  Th£N  X;r*; 

",4^? 

IP  xa  *10  THEN  PRINT  "0  = 

1 . LB'  • xa.  '  " 

~4  “5 

IP  x  a ’ 1  AND  xa  10  AND  x;n*i  " 

'*•£».  pc  ;  dp  -  :  ,  D5  .  - 

“48? 

IP  xa  1  AND  X ; r*C  THEN  OPIN' 

'CP  -2,-l..S"xa 

“4  55 

IP  ‘a*l  ThEN  PPIN"  "CP  -1,-1. 

.£1.0" 

“45? 

NE ■ *  ■  a 

"455 

PPIN"  =.  =a  ?,?" 

-  t  V 

i-*.*:  "cc ;  p  • 

z  ?  z 

' -c  •?=■'. r  ' 2  • "a •  5"£p  •«:  ec 

“*  *"  T 

Z  *  V 

•  ■  *  a-'*-»«-  *5', 

"5 : 5 

PR  I»;*  ’Pa  0 .  ,  •  - 

IF  lr»*0  ANl  Ya  0  ThEN  P=INT 

"PP  2,0" 

“525 

IF  va  1  AND  >a  0  ThEN  PRINT 

•CP  -4,-.:s.lB0.pp 

"SIC 

IF  va*0  THE*.  PPIN'  'DP  -D.-.2 

=. _E0" 

’’Sir 

IP  v  a  1  AND  Ihf.I  THEN  POINT 

•2°  -5, -.25. LB '  .Ya.  "" 

"54? 

IF  x a  - 9  AND  Ic.t-2  THEN  PRINT 

"CP  -4.-.25.LB" . Ya.  " " 

"545 

IF  Ya-0  AND  T a '  10  AND  Iht-D  THEN  PRINT  "CP  25, 

"55? 

IF  v a* 1  THEN  PRINT  CP  -4.-.D 

5  .  lS 1 . 0 " 

7555 

NEx'  Ya 

756? 

I c  Ox  d*0  then 

75S5 

XlabelS*  x - 

7570 

y labels* ” Y ' 

-■c-c 

ElSE 

7560 

BEEP 

756S 

INPvJT  "'ENTER  X  -LABEL  '  .XlabelS 

759? 

BEEP 

7595 

INPUT  "ENTER  Y-lABEl" .YlabelS 

7600 

END  IP 

7605 

PRINT  "Sc  1.5,2  Pu  PA  50,-10  CP", 

-LEN( XlabelS )/2 . "0.LB* 

.  x  e 
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0,1.  i_B " . Ylaoeit . 


'SIC 

-PINT  “PA 

-11,50  CF 

C  ,  ’  .  -LEN  Y 1 ase  1 i  )  S  »5 

’£ •  "DI  0, 1 . lB’ ; Vlaoeif 

”S  1 B 

PRINT  "CP 

0,0  DI" 

"'BSC 

Nr  jri“0 

-ESS 

BEEP 

'53C 

INFUT  " UANT  TO  PlOt 

DPT  A  FROM  A  RIlE  .  l«v 

,  £  *  N  ;  “  "  ,  G  k 

"S35 

G*  p“C 

7640 

IF  C*  *  1  T HEN 

1 4  ^ 

BEEP 

""SBC 

INPUT 

"ENTER  THE 

NAME  OF  THE  PlOT  DATA 

FILE" ,C_f : le5 

7SS5 

ASSIGN 

(JPile  TO 

■  CSC 

i.r  1COWD'  .  5  rttN 

=  sx. 

7£EE 

S>«0 

"E"0 

S,»P 

EaS-0 

"ESC 

Ex  y *0 

-0  55 

He*’ 

"S90 

BEEP 

~  =  9S 

INPuT 

“ENTER  THE 

BEGINNING  Run  number 

np  Z  s'-  M  r* 

-'■'00 

Npa :  .-s 

-9 

""PE 

EEE  = 

.  X. 

I  NP-JT 

” E N T E R  the 

NUMBER  0 -  x -  *  —  A I R S  S 

'  i  i‘ ,Nce:r: 

:  s 

BEE- 

INPUT  "SELECT  T JBE  NUMBER  ■ SsTCR . 1-SE20N2 . '.I  tube 

Nrur-  =  Nru,'  +  Nc3i''5 

pr:nte=  is  i 


~"40 

PR  In*  USING 

4  x ,  select  a  Symfc c. 

«  4  2 

PPIN"  USING 

4X  ,  "  " 

Star  2  Plus  siyr . 

'  '  -  C 

PR  IN’  USING  ■ 

4*  ,  -  " 

Circle  4  Scaa-e' 

-T  — 

CR I NT  USING  “ 

4X  ,  ■  “ 

Ronaus'"" 

"*  T  21  "* 

PRIN'  USING  ' 

4X  ,  '  ' 

R i ght -s i de-up  triangle' 

7"cE 

PRINT  USING  " 

4  X  , 

Up-s  i  ae-dowr  triangle''' 

I N- _  ~  £  y  n 

77“*- 

PAINTER  IS  70S 

"“EC 

IP  S , 1  Then 

DRIN* 

SM  • 

7  7  h  - 

Ir  s y h* s  Then 

°R  I  NT 

“  SM+  ” 

3  v 

I-  Syi— 2  Then 

PRIN* 

"SMo" 

”''?r 

M  a  =  1 1  ube*? 

”  -  i  ■> 

S  1 1 

I-  m;  i  *-,£>. 

“7  Z 
tjii 

FCR  1*0  TC 

r  I  il 

EN'ER 

@p:  Is 

7  z  •  z 

N£  X  ~  I 

£•-:  ;~ 

c 0-  I  =  i  TO  Npa : r * 


*C  _  0 

ENTEP  : I e . xa  ,  '*  a 

7SS5 

S**3»*X£ 

7340 

£  y  ■  ?  y  ♦  » a 

734E 

S*2*Sx  2*Xa"2 

7S5C 

Sxy”S>.  y  +  Xa*Ya 

-J 

CD 

in 

tn 

X * < Xa-Xrin  )  •  S  ?  * 

78E0 

Y“  '  Ya-Ypnin  >  •  S  f  y 

7£53 

IF  V-100  OR  Y<0  THEN  7910 

7S70 

IR  Sym>3  THEN  PRINT  "SM" 

7S"5 

IP  SyM<4  THEN  PRINT  "SR  1.4, 2. 4 

7080 

PRINT  "PA" tX , Y, "■ 

7885 

IP  Syf-t>3  THEN  PRINT  "SR  1.2. 1.6 

7590  IP  Sym-4  THEN  PRINT  "UC2 ,4, 99. 0,-8, -4, 0,0, 8, 4,0. " 
7S9S  IF  Syn-S  THEN  PRINT  ’ UC3 , 0 , 99 . -3 , -6 . -3 , B . 3 . 6 , 3 , -6 , 
7900  IF  Sym-6  THEN  PRINT  ” UC0 . 5 . 3 , 99 . 3 . -8 . -6 . 0 . 3 . 8 . " 
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"905 

"910 

"915 

?s:e 

7925 

7930 

7535 

•^940 

"945 

7950 

7955 

7960 

7365 

7970 

"975 

"990 

7985 

7995 

8000 

8005 

5010 

6015 

5030 

8C35 

8030 

8035 

8040 

8045 

8050 

8055 

8060 

8065 

8070 

8075 

8080 

8085 

8050 

8095 

8100 

8105 

6110 

6115 

8120 

6125 

8130 

8135 

6140 
814  = 
8150 
8155 
8160 
6165 
8170 
6175 

eise 

8185 

8190 

8195 

8200 

8205 


IF  Syi"-7  Then  PRINT  "UC0.-5.3,99,-3.8,6.e.-3,-S.  •' 
NEXT  I 
BEEF 

INPUT  “WANT  TO  LABEL  (  1 - Y , 0*N >? " . II b 1 
IF  Ilbl-1  THEN 

IF  5ynx3  THEN  PRINT  "SM" 

IF  Syr^4  THEN  PRINT  "SR  1.4,2.4- 
PRINT  "PA" ,Xal . Yal , " " 

IF  Sym/3  THEN  PRINT  "SF  1.2. 1.6“ 

IF  Sym-4  THEN  PRINT  ■ UC2 , 4 , 99 , 0 , -8 . -4 , 0 , 0 , 8 , 4 , 0 ; ■ 

IF  Syn»5  THEN  PRINT  "UC3 , 0 . 99 , -3 . -6 , -3 , 6 , 3 , 6 , 3  , -6 . * 

IF  Syn-6  THEN  PRINT  "UC0 ,5. 3, 99,3, -8, -6, 0,3, 8;” 

IF  Sym«7  THEN  PRINT  "UC0 . -5 . 3 , 99 , -3 , 6 . 6 , 0 , -2 , -8 . " 
PRINT  "SM" 

IF  Syn<4  THEN  PRINT  "PR  2,0" 

PRINT  “PP  2, -1.0. SR  1.0,1.8,LB".D_file*;-- 
Ya 1 =Y  a  1 -5 

INPUT  " UANT  TO  ADD  ANOTHER  STSINS  i 1  * Y , 0=N >~ " , I  as 
IF  Ias*l  THEN 
BEEP 

INPUT  "ENTER  THE  ETR I  No '  , Labe  1 S 
PRINT  "FP  2.0. SR  1. 0.1. 8. LB" .Labels. 

GOTO  7990 
END  IF 
END  IF 
BEEF 

INPUT  "WANT  TO  COMBINE  ANOTHER  cI_£t  (1*Y,0»N  ,  Icont 

ASSIGN  @Fj ie  TO  • 

IF  Iconc.  '0  THEN  7845 

Ils-1 

BEEP 

input  "WANT  t;  fwot  A  lEAST-SDjAPES  line  ( :-DEF«YE£,0*NC 
IF  1 1  a- 1  then 
BEEF 

INPUT  "SELECT  CURVE  TYPE  < 0-SCLIC , 1 -DASHED  1 1 1 
:  1 1 » 1 1 1  ♦  i 
PRINT  "SM ' 

IF  I e  'P'0  Th£N 

B b  =  ■ N- jn ♦ 5 - y - 5 y • 5 >  ' •  'Nrun»5'2-S- ' 2 
ELSE 

£  t  *  .  ~  S 
END  IF 

Aa=<Sy-6b*S'  l/Nrur 
PFINTEP  IS  1 

PRINT  USING  1 0 x , " " a  -  " "  ,  MZ  . 3DE "  . Aa 
PRINT  USING  "10/ .""n  -  " "  ,  M2 . 3DE  "  .  BD 
PRINTER  IS  705 
I  ^"0 

FOR  Xa*Xi*un  TO  Xna»  STEC  (Xma'-Xmin) 

Ya*Aa*xe*8b 

Y  -  <  Ya-Yrun  > »  S  f  y 

X»(Xa-Xmn  j  •  S  f 

IF  Y<0  THEN  Y-0 

IF  Y > 1 00  THEN  GOTO  8220 

IF  Ilt-1  THEN 

PRINT  "PA" , X,Y, "PD" 

ELSE 

In-In+1 
Ir»In  MOO  I  1 1 

IF  Ir* 1  THEN  PRINT  "PA" , X , Y , 'PD " 
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rj  u'j  fj  in  o  in  e*»  in  in 
.  ,  m  n  n  n  n  ^  u»  in 


IF  Ir-0  THEN  PRINT 
END  IF 
NEXT  Xa 
PRINT  •PU" 

END  Ic 
I  comb“0 
GOTO  7620 
END  IF 

PRINT  "PU  PA  0,0" 

PRINT  "PU  PA  0,0  SP0" 

6260  SUBEND 


"PA" ,X,V, "PU" 


t 
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